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ABSTRACT 


The  present  work  is  a  detailed,  theoretical  study  of  characteristic 
ignition  behaviour  of  solids,  which  are  capable  of  self  sustained  combustion, 
and  ignite  solely  due  to  gas  phase  processes.  The  physical  and  mathematical 
situation  in  which  ignition  occurs  is  that  of  a  solid  exposed  to  a  "tailored'' 
reflexion  of  a  shock  wave  from  the  proepllant  surface.  The  primary  object  of 
the  investigations  of  Part  I  is  to  analyse  the  structure  of  the  gas  phase 
adjacent  to  the  surface  of  a  homogeneous  solid  propellant  during  ignition, 
and  the  behaviour  of  the  gas  phase  throughout  the  entire  ignition  process  tc 
the  establishment  of  steady  state  burning.  It  has  been  found  —  through 
numerical  computations  —  that  a  long  time  after  the  propellant  surface  is 
exposed  to  the  ignition  stimulus,  there  exist  two  distinct  reaction  zones  in 
the  gas  phase.  An  approximate  theoretical  criterion,  chosen  to  represent 
some  gas  phase  property  within  the  secondary  diffusion  reaction  zone,  exhibits 
an  approximately  inverse  square  dependence  of  ignition  times  upon  the  external 
oxidizer  mole  fraction  as  has  been  observed  experimentally.  Asymptotic 
analyses  have  been  carried  out,  where-ever  possible,  to  check  the  results  of 
numerical  solutions.  The  agreement  was  found  to  be  excellent. 

To  facilitate  comparison  of  the  theoretical  results  with  experiments, 
and  to  explain  the  variations  present  in  the  experimental  data,  a  new  theoretical 
ignition  criterion,  based  upon  the  time  intensity  of  light  emission,  was 
developed.  The  results  obtained  with  this  criterion  are  in  excellent  agree¬ 
ment  with  those  of  experiments.  It  has  been  found  that  a  low  emission  criterion 
corresponds  to  the  development  of  a  small  temperature  rise  in  the  system  — 
dependent  upon  the  development  of  the  primary  reaction  zone  —  and  a  high 
emission  criterion  corresponds  to  a  long  time  criterion  dependent  upon  the 
development  of  the  secondary  reaction  zone.  With  the  help  of  the  new  criterion, 
the  type  of  data  needed  for  experimental  determiantion  of  the  ignition 
mechanism  of  a  given  solid  has  been  described.  From  numerical  computations, 
as  well  as  from  asymptotic  solutions,  the  marked  influence  of  the  total  pressure 
of  the  gas  phase  upon  the  ignition  behaviour  of  the  solid  has  been  described. 

Since  present  day  propellants  may  produce  a  highly  reactive  oxidizer 
upon  decomposition,  the  case  of  a  propellant  producing  an  oxidizer  much  more 
reactive  than  the  ambient  oxidizer  has  been  considered  in  Part  II.  The  results 
show  that  the  influence  of  the  ambient  oxidizer  upon  the  ignition  behaviour 
progressively  decreases  as  the  reactivity  of  the  evolved  oxidizer  increases 
leading  to  the  conclusion  that  either  very  little  oxidizer  is  decomposed 
during  the  shock  tube  ignition  of  a  typical  solid  propellant  or  that  the 
kinetic  properties  of  the  evolved  oxidizer  are  at  best  the  same  as  those  of 
the  ambient  oxidizer . 

The  ignition  of  a  composite  or  heterogeneous  propellant  is  dealt  with 
in  Part  III.  The  model  considers  two  dimensional  diffusion  and  conduction  in 
the  gas  phase  and  a  solid  consisting  of  an  annulus  of  fuel  surrounding  a 
cylinder  of  solid  oxidizer.  Here  also,  the  results  obtained  with  short  time 
temperature  criterion  and  the  emission  criterion  are  in  very  good  agreement 
with  experiments.  In  particular,  an  approximately  inverse  square  dependence  of 
ignition  times  upon  the  oxidizer  mole  fraction  with  almost  any  realistic 
ignition  criterion  is  found.  The  theory  also  indicates  that  the  ignition 
behaviuur  of  a  heterogeneous  solid  propellant  approached  that  of  a  homogeneous 


propellant  as  tha  particle  size  is  decreased. 

Thus  the  present  theoretical  studies  of  decomposable  solids  leads 
to  the  conclusion  that  a  purely  gas  phase  mechanism  alone  can  explain  the 
experimentally  observed  ignition  behaviour  of  both  homogeneous  and  hetero¬ 
geneous  propellants  under  shock  tube  conditions. 

In  Part  IV  a  brief  discussion  of  radiant  ignition  is  presented. 
Ignition  of  propellants  by  radiant  energy  has  been  widely  employed  to 
screen  the  propellants  according  to  their  ignitability.  Gas  phase  model 
describing  the  ignition  of  propellants,  including  the  effect  of  absorption 
of  radiant  energy  by  gret  gases  is  considered.  Numerical  results  agree 
with  the  approximate  analysis  of  Waldman  and  Summerfield.  However,  the 
present  treatment  does  not  conclusively  indicate  the  dominance  of  gas  phase 
mechanism  over  heterogeneous  surface  reactions  during  ignition  transients. 


iv 


NOMENCLATURE 


A  QR/nCpEbjdimensionsless  heat  reLease  parameter 

B  pre-exponential  factor  of  pyrolysis  rate 

b  dimensionless  mass  flux  parameter 

c  specific  heat  of  the  gas  phase,  O.J  cal/g°C 

2 

local  mass  diffusivity,  cm  /sec 

2  2  2 
D  =(p*0/P  '  ),  effective  mass  diffusivity,  cm  /sec. 

8 

E^,E2>E  activation  energy  of  gas  phase  reactions,  cal/mole 

E  activation  energy  of  pyrolysis,  cal/mole 

P 

®ol’®l’®2  stac^st^ca^  wef8^t  of  ground  and  excited  states,  cals 
h^  heat  of  pyrolysis,  cal/g 

I  intensity  of  electromagnetic  radiation  cm/sec 

J  relative  response 

E  R 

M  u  ,  a  dimensionless  quantity 

C  kE 

,  2 

m  total  mass  flux  issuing  from  the  propellant  surface,  g/cm  -sec 

m^  initial  mass  flux,  g/cm^-sec  (at  400°K) 

n^,n2>n  molar  stoichiometric  ratio  of  gas  phase  reaction 

Q1>Q2>Q  heat  of  combustion,  cal/mole 

R  universal  gas  constant,  1.986  cal/mole  °K 

r  radial  distance 

t  time,  sec 

ic 

t  ignition  delay,  millisec. 

o 

T,Tq  local  initial  gas  phase  temperature,  K 

T  ,T  ,  local  initial  solid  phase  temperature,  °K 
s  s  o  ^ 

V  local  velocity,  cm/ sec 


linear  distance,  cm 


<a 


V 


yf>ybx 


y"  >y 

ox  ox 


2,  ,Z0,Z 


€1,€2 


w%» 


K  >  X 

S  g 


P’PgPs 


P2>pi 


respectively,  mole  fractions  of  fuel  and  evolved  oxidizer 
gas  phase  supplied  and  total  oxidizer  mole  fractions 

(y  =  y'  +  y"  ) 

'ox  'ox  'ox 

initial  oxidizer  mole  fraction  in  the  gas  phase 

mole  fraction  of  the  inert  gases 

pre-exponential  factor  of  gas  phase  reactions 

mf  ,  ratio  of  fuel  mass  flux  to  total  mass  flux 
m 

2 

thermal  diffusivity  of  solid  phase,  cm  /sec 

2 

thermal  diffusivity  of  gas  phase,  cm  /sec 

_ v  ,  dimensionless  heat  of  pyrolysis 

c  E 
P 

ignition  criteria 
energies  of  the  excited  states 
dimensionless  reactant  concentration 
inert  gas  concentration,  dimensionless 
dimensionless  gas  phase  temperature 
dimensionless  solid  phase  temperature 
initial  gasphase  temperature,  dimensionless 
s( ,  a  dimensionless  quantity 
stability  parameter 

thermal  conductivities  of  solid  and  gas  phase 
ratio  of  initial  gas  to  solid  phase  densities 
dimensionless  distance  in  gas  phase 
dimensionless  distance  in  solid 

local  gas  phase,  initial  gas  phase  and  solid  phase  densities 
dimensionless  radial  distance 
dimensionless  outer  and  inner  radii 


ratio  of  solid  to  gas  phase  thermal  dif fusivities 
dimensionless  time 

extinction  coefficient  for  radiant  transmissions  cm 
Boltzman  constant 

solid 
oxidizer 
condensed  phase 
gas  phase 
initial 
propellant 
at  infinity 


ignition  condition 
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INTRODUCTION 


Ignition  of  solid  propellants  has  been  a  topic  of  considerable 
interest  in  the  development  of  operational  motor  systems  for  some  time.  In 
critical  applications,  such  as  missiles  wh^re  a  large  amount  of  energy  has 
to  be  released  in  a  very  short  time,  it  is  obvious  that  the  ignition  system 
has  to  operate  repeatably  and  reliably.  Design  of  such  igniter  system  is 
eased  by  some  knowledge  of  ignition  mechanism. 

Conceivably,  to  a  certain  degree,  the  ignition  mechanism  depends 
upon  the  ignition  system  employed.  For  example,  if  the  ignition  is  achieved 
by  the  application  of  a  hypergol,  like  chlorine  trifluoride,  the  run-away 
reactions  occur  at  the  propellant  surface  and  are  termed  as  hetrogeneous  re¬ 
actions.  There  are  several  methods  of  igniting  a  solid  propellant  but  the 
mechanism  of  ignition  will  fall  into  one  of  the  following  categories  (1) 
solid  phase  reactions,  (2)  hetrogeneous  reactions  and  (3)  gas  phase  re¬ 
actions.  Determination  of  the  exact  mechanism  by  which  a  propellant  ignities 
is  extremely  complex  in  any  case. 

However,  in  certain  cases,  like  ignition  in  shock  tubes,  it  is 
argued  that  ignition  occurs  by  gas  phase  reactions.  This  is  understandable 
because,  the  surface  temperature  reached  in  shock  tube  experiments  are  too 
low  for  any  significant  surface  reactions  to  occur.  In  the  case  of  hyper- 
golic  ignition,  the  run-away  reactions  are  indisputably  hetrogeneous.  But 
in  arc  image  ignition  situations,  the  mechanism  by  which  the  ignition  occurs 
is  not  obvious.  In  steady  burning,  solid  phase  reactions  appear  to  be 
unimportant  in  the  case  of  composite  propellants.  But  it  is  not  clear  whether 
the  ignition  is  due  to  hetrogeneous  reactions  or  homogeneous  gas  phase  reactions. 
Possibly,  it  is  a  combination  of  both.  Theoretical  models  employing  purely 
gas  phase  or  purely  heterogeneous  reactions  predict  nearly  the  same  behaviour, 
and  no  undisputably  deterministic  experiments  have  been  performed.  In  the 
absence  of  proof  for  any  particular  mechanism  controlling  ignition,  the  only 
way  out  for  an  investigator  is  to  examine  theoretically  the  ignition  behaviour 
assuming  that  the  propellant  ignition  occurs  by  one  of  the  above  mentioned 
mechanisms . 

The  object  of  the  present  analysis  is  to  establish  the  ignition 
behaviour  based  on  the  assumption  that  the  run-away,  thermo-chemical  reactions 
leading  to  ignition,  occur  in  the  gas  phase  adjacent  to  the  propellant  surface. 

The  beginning  of  the  history  of  ignition  theories  for  solid  pro¬ 
pellants  was  marked  by  a  simple  solid  phase  ignition  theory  postulated  by 
Hicks.  In  his  theoretical  model  Hicks  assumed  that  ignition  was  caused  by  run¬ 
away  thermo-chemical  reactions  occurring  in  a  thin  region  of  the  solid  adjacent 
to  the  propellant  burning  surface.  The  entire  ignition  lag  was  composed  or 
heat-up  period;  the  time  required  to  establish  a  gas  phase  flame  was  assumed  small. 
The  theoretical  predictions  of  the  ignition  behaviour  were  in  good  agreement 
with  the  experimentally  observed  behaviour  of  the  double  base  propellants 
employed  at  that  time.  However,  it  has  been  found,  that  in  the  case  of  modern 
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day  composite  propellants,  the  ignition  behaviour  depends  considerably  upon 
the  properties  of  the  adjacent  gas  phase;  solid  phase  reactions  are  not 
important  during  rapid  ignitions.  Conclusive  proof  of  this  statement  came  ^ 
from  shock  tube  experiments  on  both  double  base  and  composite  solid  propellants  . 
In  addition,  experiments  on  steady  burning  of  composite  propellants  have,  in 
the  past,  indicated  that  vitally  important  heat  release  occurs  in  the  gas  phase. 
These  and  other  considerations  led  to  the  study  of  models  involving  the 
participation  of  gas  phase  processes  during  ignition  transients. 

Studies  made  at  Princeton  University  during  the  period  from  1958  to 
1963,  on  several  double  base  and  composite  propellants,  marks  the  beginning 
of  intensive  research  in  solid  propella.it  ignition.  Based  on  certain  experi¬ 
mental  observations,  McAlevy  and  Summerfield  developed  a  gas  phase  ignition 
model  in  which  a  hot  oxidizing  gas  reacts  with  fuel  vapours  issuing  from  the 
propellant  surface  due  to  decomposition  of  the  solid.  Importance  was  given  to 
a  special  case  of  an  environment  containing  an  oxidizer,  as  the  propellants 
tested  would  not  ignite  unless  the  environmental  gas  contained  sufficient 
oxygen  under  the  experimental  conditions  used.  No  analytical  study  encompassing 
a  more  general  case  of  propellant  decomposing  into  fuel  and  oxidizer  vapours 
was  made.  Initial  theoretical  calculations  using  gas  phase  ignition  theory  for 
solid  fueld  showed  qualitative  agreement  with  the  observed  dependence  of 
ignition  of  composite  propellants  upon  the  environmental  gas. 

Intensive  research  activity  in  the  field  of  solid  propellant  ignition 
is  perhaps  better  illustrated  by  a  brief  mention  on  the  emergence  of  hypergolic, 
hetrogeneous  ignition  theory  developed  by  Anderson  and  Brown  2  almost 
simultaneously  to  the  development  of  gas  phase  theories.  This  theory  was 
originally  evolved  from  the  studies  on  ignition  of  solids  by  chemical  activation 
of  the  propellant  surface  by  powerful  oxidizers  like  fluorine  and  chlorine  tri- 
fluoride.  Accordingly,  it  was  assumed  in  the  hypergolic  ignition  model  that 
ignition  was  initiated  by  heterogeneous  surface  reactions  between  the  hypergol 
and  the  fuel  surface.  It  was  argued  by  Anderson  and  Brown  that  ignition  under 
shock  tube  conditions  might  be  due  to  surface  reactions  similar  to  hypergolic 
ignition  reactions,  since  theory  predicted  a  strong  dependence  of  ignition 
behaviour  upon  the  properties  of  the  adjacent  gas  phase  as  was  found  from 
experiments . 

Though  the  predictions  of  hypergolic  model  are  in  good  agreement  with 
experiments,  it  is  somewhat  dubious  that  surface  reactions  may  have  the  strong 
hold  in  initiating  ignition;  the  reported  shock-tube  generated,  propellant 
surface  temperatures  were  low,  and  oxygen  is  not  active  enough  to  cause  vigorous 
surface  reactions  at  these  temperatures.  Therefore  further  detailed  study  of 
gas  phase  processes  of  solid  propellants  (not  pure  fuels)  —  discussed  in  the 
following  sections  —  was  necessary  before  making  further  conclusions  regarding 
the  dominant  ignition  processes  of  solid  propellants.  The  present  investigation 
is  an  effort  directed  towards  better  understanding  of  the  gas  phase  processes. 
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THEORY  OF  GAS  PHASE  IGNITION  OF  HOMOGENEOUS 
AND  HETEROGENEOUS  SOLID  PROPELLANTS 


CHAPTER  I 


1 . 1  Review  of  Previous  Work: 

Theoretical  studies  on  ignition  can  be  traced  back  to  the  early  works 
of  Semenov  and  Frank-Kamenetski  who  considered  ignition  as  a  thermal  run-away 
in  a  system  consisting  of  combustible,  homogeneous  mixture  of  a  fuel  and  an 
oxidizer  in  the  vapour  phase.  In  this  type  of  ignition  it  is  assumed  that 
exothermic  chemical  reactions  are  occurring  everywhere  in  the  system  simul¬ 
taneously;  the  gaseous  mixture  ignites  after  a  short  induction  period  under¬ 
going  a  small  temperature  rise.  A  system  in  which  the  reactants  are  premixed 
is  called  a  homogeneous  system.  Contrasted  to  the  ignition  of  a  homogeneous 
system  is  the  ignition  of  heterogeneous  systems  in  which  the  reactants  are 
initially  unmixed;  the  reactants  may  exist  in  different  phases.  A  detailed 
discussion  of  homogeneous  and  hetergeneous  reactions  is  given  in  reference  (3). 
However,  a  few  basic  differences  between  homogeneous  and  heterogeneous 
ignitions  may  be  worth  mentioning:  (1)  thermal  ignition  of  a  homogeneous  system 
depends  entirely  upon  the  thermo-chemical  properties  of  the  mixture  whereas  the 
ignition  of  a  heterogeneous  system  depends  upon  the  rate  of  consumption  or  the 
rate  of  supply  of  one  or  either  of  the  reactants;  (2)  theoretical  determination 
of  the  ignition  of  a  homogeneous  system  does  not  depend  upon  the  choice  of 
ignition  criteria  whereas  that  of  a  heterogeneous  system  very  strongly  depends 
upon  the  choice  of  ignition  criteria  as  discussed  in  later  chapters.  In  what 
follows  only  heterogeneous  systems  in  which  both  reactants  exist  in  the  same 
(gaseous)  phase  are  considered. 

The  first  attempt  to  describe  mathematically  and  physically  the 
ignition  of  heterogeneous  systems  was  made  by  McAlevy  and  Summerfield.  In 
their  simplified  version  of  the  gas  phase  model,  they  assumed  that  the  fuel 
vapours  after  decomposition  at  the  proepllant  surface  diffused  outwardly  from 
the  surface  into  the  adjacent  gas  phase  containing  a  hot  oxidizer.  Vigorous 
chemical  reactions  occur  in  the  gas  phase  with  the  release  of  large  amounts  of 
heat  leading  to  thermal  ignition.  Several  simplifying  assumptions  were  made  in 
their  theoretical  analysis  to  solve  the  mathematical  ignition  model.  The 
ignition  event,  according  to  McAlevy,  was  defined  as  the  instant  at  which  the 
heat  generation  at  any  location  in  the  gas  phase  exceeds  the  heat  losses  at 
that  location.  Defining  ignition  in  this  manner  was  convenient  for  mathematical 
analysis,  as  well  as  being  physically  intelligible.  The  model  was  not  an 
unequivocal  success,  however,  though  the  predictions  were  in  rough  agreement  with 
the  observed  behaviour  of  solid  propellants  under  shock  tube  experiments.  The 
theory  indicated  that  t  a.  (y qX'  ^  whereas  the  experiments  showed  a  much 
stronger  dependence  of  ifnition  times  upon  the  initial  oxidizer  mole  fraction. 
But,  before  discarding  the  gas  phase  theory  of  ignition,  it  was  necessary  to 
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remove  some  of  the  restrictive  assumptions  of  the  McAlevy/Summerf ield  model. 

Later  work  by  Hermance3  included  a  more  complete  model  in  which  simultaneous, 
transient,  mass  diffusion  and  chemical  reactions  occurred  within  the  gas  phase; 
both  the  solid  surface  temperature  and  the  surface  temperature  dependent 
pyrolysis  rates  were  assumed  to  be  constant  in  time.  The  model  demonstrated 
the  strong  effect  of  neglecting  the  consumption  of  the  reactants  (McAlevy's 
model  neglected  reactant  consumption)  and  the  dependence  of  ignition  delay 
on  the  ignition  criteria. 

The  mathematical  criterion  for  defining  the  ignition  event  in 
Hermance's  work  was  the  first  attainment  of  a  maximum  gas  phase  temperature  which 
was  50%  in  excess  of  the  initial  gas  phase  temperature.  It  was  found  that 
McAlevy's  criterion  was  inappropriate.  Indeed  it  was  shown  by  Hermance^  that 
McAlevy's  criterion  was  satisfied  much  too  soon  in  the  transient  gas  phase 
processes  leading  to  ignition.  The  model  gave  results  in  better  agreement  with 
experiments  on  the  ignition  of  pure  fuels,  but  still  could  not  explain  the 
strong  dependence  of  ignition  delay  upon  the  initial  oxidizer  mole  fractions  for 
propellants.  The  assumption  of  constant  wall  conditions  was  suspected  to  be 
responsible  for  this  lack.  Accordingly,  self  determining  surface  conditions 
(termed  'feed  back'  conditions)  were  introduced  by  Hermance^  into  the  basic 
ignition  model.  This  refinement  established  at  least  one  interesting  aspect: 
that  even  for  fuels  having  low  volatility,  or  low  pyrolysis  rates,  at  the 
interface  temperatures  initially  established,  the  ignition  occurred  in  a  very 
short  time  (order  of  milliseconds)  which  was  in  agreement  with  experiments  on 
some  less  volatile  (e.g.  epoxy)  fuels.  Surface  temperature  dependent  pyrolysis 
rate  of  the  form  m  =  PsBe-E/RTS  was  assumed.  The  ignition  criteria  employed 
remained  the  same  as  before.  Since  the  mathematical  model  included  the  ignition 
of  pure  fuels  only,  the  work  was  considered  incomplete.  Further  improvement  — 
the  case  of  a  solid  which  after  decomposition  produces  both  fuel  and  oxidizer 
vapours  —  was  imminent. 

A  gas  phase  model  in  which  the  propellant  decomposed  into  fuel  and 
oxidizer  vapours  was  considered  by  Hermance  and  Kumar^>5,  The  model  revealed 
two  important  things:  (1)  the  propellant  ignities  in  a  neutral  atmosphere 
utilizing  its  own  oxidizer  and  (2)  nearly  steady  bruning  (or  quasi-steady) 
conditions  were  established  in  the  gas  phase  after  a  long  time  (order  of  few 
milliseconds).  The  model  also  indicated  that  t^gn  ^  (yox^_T1  w^ere  n  £  2/ 3 
depending  upon  the  composition  of  the  propellant  (a  high  percentage  of  oxidizer 
in  the  solid  showed  a  weaker  dependence  of  ignition  on  the  external  oxidizer). 

The  ignition  criterion  employed  was  the  same  as  the  one  in  references  3  and  4. 

There  were  several  obvious  areas  in  which  this  model  could  be 
improved.  These  are  (1)  incorporating  natural  boundary  conditions  -  i.e.  - 
feedback,  (2)  development  of  an  ignition  criterion  which  closely  corresponds 
to  the  experimental  ignition  criterion  of  light  emission,  (3)  use  of 
temperature  dependent  gas  phase  physical  properties  and  (4)  the  effect  of 
propellant  inhomogeneuity.,  A  detailed  theoretical  analysis  of  a  gas 

phase  ignition  model  with  the  introduction  of  the  aforementioned  aspects  is 
presented  in  the  succeeding  chapters  —  with  interesting  and  significant 
results . 


IGNITION  OF  HOMOGENEOUS  SOLID  PROPELLANTS  UNDER 
SHOCK  TUBE  CONDITIONS 


CHAPTER  II 


2.1  Physical  Model: 

The  model  considers  one  dimensional  diffusion  and  simultaneous,  gas 
phase  chemical  reaction  of  the  gaseous  species  produced  by  the  thermal  de¬ 
composition  of  a  homogeneous  propellant  under  shock  tube  conditions.  A 
schematic  representation  of  the  gas  phase  ignition  model  is  shown  in  Figure  1. 
Ignition  in  a  shock  tube  can  be  described  as  follows:  A  high  intensity  shock 
wave  is  made  to  propagate  through  a  homogeneous  mixture  of  an  oxidizer  and  a 
neutral  gas  (generally  nitrogen)  in  desired  proportions,  by  the  rupture  of  a 
diaphragm  as  described  in  reference  1.  Upon  reflection  from  the  propellant 
surface  (located  flush  with  the  end  wall  of  the  shock  tube)  the  shock  wave 
creates  a  high  temperature  and  high  pressure,  stagnant  gaseous  atmosphere  in 
the  gas  phase.  Simultaneously,  the  propellant  surface  temperature  jumps  to  a 
certain  value  determined  by  the  thermal  properties  of  the  gas  and  the  solid. 
The  propellant  surface  starts  to  decompose  at  a  rate  which  is  dependent  upon 
the  mixture  of  fuel  and  oxidizer  vapours,  which  in  turn  diffuse  into,  and 
chemically  react  with,  the  external  oxidizer  initially  present  in  the  gas 
phase.  A  major  portion  of  the  heat  released  by  this  chemical  reaction  is 
used  up  in  raising  the  temperature  of  the  gas  phase;  the  remaining  heat  is 
fed  back  to  the  solid  which  in  turn  raises  the  surface  temperature,  providing 
a  boot-strap  effect.  The  following  assumptions  were  made  in  formulating  the 
present  gas  phase  model: 


1)  The  molecular  weights  of  gas  and  solid  are  constant  and  equal. 


2)  The  propellant  is  a  homogeneous  mixture  of  fuel  and  oxidizer 

and  decompose  according  to  the  pyrolysis  law  m  =  p  B  exp(-E  /RT  )  * 

producing  a  homogeneous  gaseous  mixture  of  fuel  an§  oxidize?. 

3)  The  specific  heat  and  thermal  conductivities  of  gas  and  solid 
are  constant  and  equal. 
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4)  The  Lewis  number  is  unity  and  p  is  a  constant  independent 
of  temperature. 

5)  pZ  is  a  constant  independent  of  temperature. 


6)  The  chemical  reaction  between  the  reacting  species  is  of  global 
second  order  of  the  type 


chem 


I  ._£*)  =-fl 

*  8t'chem  Q 


£T  \ 
3t  / 


=  -  pZ 


chem 


^f^ox 


exp(-E/RT) 


with  a  constant  molar  stoichiometry  of  the  form 


fuel |  +  n  j oxidizer  j-*(n+l)  | products! 
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Assuming  constant  molar  stoichiometry  eliminates  the  necessity  of 
considering  Stephan  flow  in  the  gas  phase;  it  i  »  .ilso  a  necessary  consequence 
of  assumption  (1).  It  is  interesting  to  note  teat  assumptions  (1) , (3)  and  (4) 
are  commonly  employed  in  the  combustiot  literature  and  are  often  justified.  In 
so  far  as  assumption  (2)  is  concerned,  it  has  been  found  that  the  surface 
dependent  decomposition  rate  of  severaJ  nitrate  esters  follow  the  Arhenius 
pyrolysis  law.  Assuming  a  constant  p 2d  in  the  gas  phase  is  analogous  to  the 
assumption  of  constant  py  in  gas  dynamics  problems,  or  the  Chapman-Rubesin 
approximation.  In  the  present  case,  coupled  with  the  Howarth  transformation, 
it  reduces  the  problem  of  variable  density  into  one  of  constant  density.  The 
jump  temperature  calculated  using  the  modified  equations  is  the  correct  "double 
jump"  value.  The  jump  temperature  calculated  without  making  this  assumption  is 
only  one  half  of  the  experimentally  measured  valued.  Detailed  treatment  of 
assumptions  (4)  and  (5)  are  carried  out  in  reference(19)  where  it  is  argued 
that  the  variation  in  e_E'E^  is  much  greater  than  the  variation  of  T  alone  ana 
thus  the  assumption  that  pZ  is  a  constant  is  perhaps  justified.  Since  normal 
gas  phase  reactions  generally  involve  bimolecular  collisions,  assumption  (6) 
also  is  proper.  Equimolar  stoichiometry  may  be  true  when  a  polymer  decomposes 
into  light  weight  gaseous  hydrocarbons. 

The  present  analysis  is  greatly  simplified  by  assuming  single  step 
reaction  kinetics.  It  can  be  said  with  great  certainty  that  ignition  reactions 
are  not  that  simple,  and  single  step  kinetics  perhaps  is  a  gross  simplification 
of  the  real  physical  situation.  It,  however,  can  be  shown  that  some  multiple 
step  reactions  can  be  simplified  (mathematically)  to  single  step  kinetics  and 
since  in  most  of  the  cases  the  exact  number  of  reaction  steps  are  not  known, 
the  assumption  of  single  step  kinetics  appears  to  be  quite  reasonable. 
Dissociation  of  the  products  should  be  taken  into  account,  at  high  temperatures 
in  reality,  but  for  simplicity  dissociation  of  the  products  is  neglected  in 
the  present  analysis. 


2.2  Mathematical  Formulation: 

The  time  rate  of  change  of  concentration  of  any  species  within  the 
control  volume  is  equal  to  the  algebraic  sum  of  the  rates  at  which  the  species 
is  entering  the  control  volume  through  diffusion  and  convection,  and  the  rate 
at  which  the  species  is  consumed  or  liberated  due  to  chemical  reactions.  If  we 
assume  that  the  evolved  oxidizer  has  the  same  physical  and  kinetic  properties  as 
the  free  stream  oxidizer  equations  for  the  conservation  of  species  and  energy  can 
be  written  as  follows: 

Let  an  operator  K  be  defined  as: 


Ks  fE  + 


—  d  _1  9  <•  ~  3  \ 
p  3x  “  p  3x'P  3x' 


K(yf )  =  £  K(yox)  =  -^K(T)  =  -  pZyfyox  exp(-E/RT) 


Then 
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k  ( )  so 


2.2 


m 

r*  3x 


2.3(a) 


It  wa 3  assumed  in  the  above  equations  that  yQy,  =  y'  +  y"ox  where  y'ox 
refers  to  evolved  oxidizer  mole  fraction  and  yQ "  refers  to  free  stream 
oxidizer  mole  fraction. 


The  boundary  and  initial  conditions  for  equations  2.1  and  2.3  are 
written  as  follows: 


mw  f D 

*  <*x 
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2.3(b) 
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Wherever  needed,  the  mole  fraction  of  products  can  be  obtained  from 
y  *  1  -  (y^  +  yQx  +  y^)  which  indeed  satisfies  the  relationship 

SC. 

.  „  f  PC’Cjt),.. 

Introducing  Howarth  length  transformation  of  the  form 
and  assumptions  (4)  and  (5)  and  defining  dimensionless  constants  and  variables 

PC  • 


m0  •  fsBexK- 


6)  =  a 


2  =  «s/0 

A  e 

b  =  r*Zt>/rc>Z 

A  a  /nCfcii 

(3  Bit#  jcfB 

P  -  ty/ft 
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3  constants  and  variables 

©c  Rt/e 

05  a  ftT*/6p 

^ok*  b*°* 

7i  = 

6  *  H  ”“/$*) 


and  letting 


transforms  equations  2.1  to  2.4 

Hi1f)  =  «(»)„)  --  -  Jf  H(e;  '/*) 


and 


J»i  *  b“,  J&  =  £g* 
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subject  to  the  following  initial  and  boundary  conditions: 


o )e  rff(co,’c)o0 

^(0>t)  “tS 1 

X,(to)  --  %(«■*}  -  £ 

T-,  U.o)  «  «■}“ 

H 

0  (s>0)  -  6(">T)*6o 

e»(o,x) 

6s(-S,0  )  - 

A  ,xj 

r14  <Js,  ’  <4 
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Definition  of  5  and  ?  in  the  above  fashion  allows  one  to  calculate 
the  correct  "double  jump"  temperature  at  the  propellant  surface  when  the  equations 
are  integrated  numerically. 

Of  all  the  dimensionless  groups  listed  above,  the  most  important  ones 
are  'b*  which  is  the  inverse  of  the  first  Damkohler  number,  A  the  dimensionless 
heat  release  parameter,  and  a  the  propellant  composition  parameter.  The  parameter 
(b)  can  also  be  written  as  <o*£)  J  | 

rz 

which  is  the  ratio  of  characteristic  diffusion  times  to  characteristic  chemical 
reactions  time.  Thus  a  large  value  of  b  indicates  dominance  of  diffusional 
process  compared  to  chemical  processes  in  controlling  ignition  and  vice-versa 

For  a  given  propellant,  ignition  in  a  shock  tube  operating  under  given 

conditions,  the  characteristic  ignition  delay  x*  is  found  to  be  a  function  of 

y°°  alone.  Thus, 
ox  ’  —  * 

Z  </> 
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Experiments  show  a  dependence  of  x*  on  y«  as  x*  ^  |y»  |  where  1  _<  6  <_  3; 
theory  generally  gives  values  of  5  to  be°Setween  2/3  an§X2  depending  upon  the 
choice  of  ignition  criterion  and  the  completeness  of  the  theory. 


2 . 3  Ignition  Criteria 


Discussion  of  several  ignition  criteria  used  in  previous  gas  phase 
models  are  given  in  references  3  and  4.  In  all  these  cases,  ignition  was 
attributed  to  some  state  of  gas  phase  distributions  during  the  ignition  period. 
In  the  past,  three  ignition  criteria  have  been  extensively  employed  for 
defining  the  ignition  event.  In  his  gas  phase  model  McAlevy  assumed  that 
ignition  occurred  whenever  heat  generation  at  any  point  in  the  gas  phase 
exceeds  heat  losses  at  that  point.  This  is  mathematically  equivalent  to 
stating  that  1Z  =  q  at  ignition.  Hermance  argued  that  this  criterion  is 
not  appropriate  since  this  criterion  was  satisfied  much  earlier  in  the 
ignition  period,  and  need  not  necessarily  lead  to  incipient  ignition.  In 
his  calculations,  Hermance  employed  another  criterion  which  stated  that 
ignition  occurred  when  the  maximum  gas  phase  temperature  at  any  station  in 
the  gas  phase  exceeds  the  initial  gas  phase  temperature  by  a  specified  amount. 
Thus,  at  ignition,  T*  =  a  T  where  a  >  1.  For  reporting  the  numerical 
results  Hermance  used  a  value  of  a  =  1.5.  Statement  of  ignition  in  the  above 
terms  is  equivalent  to  specifying  a  certain  rate  of  temperature  rise  in  the 
gas  phase. 


Another  Ignition  criterion  often  employed  is  the  integrated  heat 
balance  criterion.  According  to  this  criterion,  ignition  is  defined  as  that 
instant  at  which  the  heat  generation  due  to  chemical  reactions,  integrated 
over  the  entire  gas  phase  at  a  given  time  is  equal  to,  and  increasing  more 
rapidly  than,  the  heat  conducted  out  of  the  boundary  layer  at  the  solid-gas 
interface.  Mathematically  expressed,  ignition  occurs  when 

Q(  generated)  +  Q  (conducted)  ^  0  t  =  t^ 

.  . 

or  0.  jVztyV*K-  VWctx  +  Aj  ~  )  >o 

o  *  & 

Theoretical  work^  indicated  that  this  criterion  gives  results  close  to  McAlevy* s 
criterion. 


All  the  three  criteria  discussed  above  appear  to  be  physically  reasonable 
and  constitute  several  approximate  mathematical  basis  for  defining  ignition. 
However,  it  should  be  noted  that  these  criteria  have  nothing  in  common  with  usual 
experimental  ignition  criteria,  or  so  it  appears.  Unless  a  relationship  can  be 
established  between  a  theoretical  and  an  experimental  ignition  criterion, 
theoretical  predictions  are  of  limited  significance. 

2 . 4  Asymptotic.  Analysis  at  Short  Times: 

Since  all  the  above  mentioned  ignition  criteria  correspond  to  ignition 
times  very  close  to  the  beginning  of  the  ignition  period,  certain  approximations 
can  be  made: 
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1.  Convective  transport  terms  can  be  neglected  compared  to 
diffusive  transport  terms  in  equations  2.6  to  2.8. 

2.  Propellant  surface  temperature  can  be  assumed  constant 
with  respect  to  time. 

3.  Reactant  consumption  up  to  ignition  can  be  neglected. 

3 

Numerical  computations  by  Hermance  ,  have  shown  that  neglecting 
convective  terms  does  not  introduce  any  significant  error  in  the  final  results. 
Although  it  can  be  shown  without  much  effort  that  assumption:  (1)  leads  to  a 
gross  violation  of  the  conservation  law,  convection  is  ignored  in  the  approximate 
analysis.  Assumption  (2)  is  well  justified  for  high  initial  pyrolysis  rates, 
since  feed  back  is  not  operative  under  these  circumstances.  Assumption  (3)  is 
also  valid  since  for  large  initial  oxidizer  mole  fractions  and  high  pyrolysis 
rates,  reactant  consumption^was  indeed  found  negligible.  With  these  it  can 
be  shown  that  4,5  T*  ^  (wAb  y“x)  £or  0Xygen  rich  and  t*  v  (nb^w  A  ctl-a) 
for  neutral  atmospheres  respectively.  These  results  are  in  good  agreement  with 
the  behaviour  of  numerical  solutions . 


2.5  Asymptotic  Analysis  at  Long  Times: 

At  long  times  (order  of  few  milliseconds) ,  quasi-steady  conditions 
prevail  in  the  gas  phase.  There  are  two  distinct  reaction  zones.  First  is 
a  premixed,  primary  reaction  zone  of  few  microns  thick  adjacent  to  propellant 
surface.  The  conditions  in  the  primary  zone  are  unchanging  and  the  zone  remains 
fixed  in  space.  There  is  a  second,  diffusion  reaction  zone,  remote  from  the 
propellant  surface  which  moves  with  a  constant  velocity  relative  to  the  pro¬ 
pellant  surface.  The  conditions  in  the  secondary  reaction  zone  are  slowly 
varying  functions  of  time.  Since  the  behaviour  of  secondary  reaction  zone 
depends  very  strongly  upon  the  external  oxidizer,  and  we  are  only  interested 
in  the  influence  of  external  oxidizer  upon  ignition,  we  only  need  to  consider 
the  secondary  reaction  zone.  This  is  facilitated  by  a  change  of  co-ordinates 
and  then  examining  the  secondary  reaction  zone  with  respect  to  this  new 
co-ordinate  system  in  which  the  reaction  zone  remains  stationary.  Assuming 
that  concentration  profiles  of  oxidizer  and  fuel  vapours  can  be  approximated  by 
power  law  expressions  ,  the  following  relationship  holds  true  in  the  secondary 
reaction  zone: 

37*  *  ,/t} 

»  W  *  -  |  >  -  -  -  y  ^  6  V  v  V  r  v  # 

-  ei-M 

X  --  ‘z  *  ">»  %y 


2.10 


Here  t  is  measured  from  the  instant  at  which  a  strong  diffusion  reaction  zone 
is  formed.  The  constants  c. ,0^,  m^  and  m^  can  be  determined  by  the  require¬ 
ments  that  at  5  =  0,  iiQ  =0;  C  =  6,  n  x~  n“x>  =  0  and  on  the  average 
stoichiometry  is  maintained  in  the  seconHary  zone. 

Ihus’  VCL'-(V</] 

’CtCV^J 

Close  to  steady  state,  the  total  rate  of  heat  generation  from  the 
secondary  reaction  zone  is  approximately  equal  to  the  rate  of  heat  loss,  which 
is  approximately  constant.  Therefore: 

v*  «€  i&h-ufw-'Ws  ,12 


letting  VW  we  get 


0 


* 


4  c  £_ 


Since  the  integral  in  the  above  equation  is  approximately  constant,  we  obtain 
6  la  1/An°°  •  Consequently  the  energy  equation  takes  the  form 

O^C 


subject  to  the  condition^  that20(O,T).  =  0.,  0($,t)  =  0q  and  0(5,0)  =f(£/6). 

If  we  define  t  and  5  as  t  =  t6x  and  5  s  ?o  to  completely  eliminate  6  from 
the  equation  2.14  and  the  corresponding  boundary  conditions,  the  solution  then 
is  of  the  form  0  =  0(0^,0  ;t,,5)  Thus  for  any  suitable  choice  of  ignition 
criteria,  the  ignition  t*  ?s  given  by: 


(m:7 


Since  n00  =  b  y®  and  A  'v  ~,  it  can  be  expected  that  at  long  times  t*  (y«°  ) 

and  t*  2X  ±  (m°|,  b  °X 

The  above  analysis  also  applies  to  pure  fuels  igniting  in  oxygen 
rich  atmospheres,  since  the  effect  of  any  ozidizer  supplied  by  the  solid  is 
obliterated  by  the  presence  of  primary  reaction  zone  where  the  oxidizer  from 
the  solid  is  consumed,  leaving  pure  fuel  vapours  to  react  with  the  external 
oxidizer.  Furthermore,  under  conditions  of  small  y m  ,  where  the  primary  zone 
is  dominant,  the  long  time  ignition  delay  is  also  independent  of  mQ. 

For  a  propellant  igniting  in  a  neutral  atmosphere,  the  power  law 

expressions  for  ancj  n 

£  ox  become 
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Furthermore,  close  to  steady  conditions  in  the  gas  phase, 

.  I  «  > 


/**  j  0-^)  e*K-  '/©> 


.  A  _  £_  2.17 

-  Qtosi  j[ 

1/2 

This  gives  the  relationship  6^6  •  Using  this  and  substituting  2^15  and 

2.16  into  energy  equation  and  redefining  t  and  ?  as  t  =  t  and  E,  =  ?6,  it 

can  be  verified  that  the  energy  equations  becomes  independent  of  all  constants, 

from  which  the  relationship  r*  ^  b-^-  and  x*  n,  ct(l~a).  Since  b  V— )2  and 
O  *  in 

r  ^  m^-,  even  for  a  propellant  igniting  in  a  neutral  atmosphere,  the  ignition 

time  is  independent  of  volitility. 

In  the  above  analysis  the  assumption  was  that  the  reaction  zone  thick¬ 
ness  depended  upon  heat  losses.  This  is  in  contrast  to  that  made  in  estimating 
the  thickness  of  a  steady  state  diffusion  flame.  In  the  former,  one  obtains 
Svo  ‘/n«  whereas  the  latter  gives  The  numerical  computations 

indicate  that  **  YvSSl  is  nearly  true;  in  effect  then,  it  can  be  concluded 
that  finite  reaction  rate  kinetics  play  an  important  role  in  the  establishment 
of  the  steady  state. 

In  the  preceding  derivation,  no  mention  was  made  of  any  specific 
ignition  criterion  for  the  theoretical  definition  of  ignition.  The  ignition 
criterion  could  be  either  Q  *  =o<  0O  ,  where  or  o.  For  presenting 

numerical  result  a  criterion  of  the  type  ‘t/0rn  0  was‘  chosen.  In 

reference  (8)  it  is  shown  that  this  also  is  an  appropriate  criterion. 

* 

The  key  assumption  in  presenting  the  long-time  analysis,  was  that 
the  emission  criterion  employed  in  experimental  detection  of  ignition  must 
in  some  sense  be  related  to  the  establishment  of  steady  conditions.  In  other 
words,  emission  of  light  during  propellant  ignition  is  delayed  until  the 
establishment  of  steady  burning  conditions.  Whether  this  is  true  is  not  known. 


2.6  Emission  Criterion 

Theoretical  studies  are  of  limited  significance  unless  they  offer 
a  means  of  comparison  with  experiments.  Though  all  the  gas  phase  theories 
discussed  previously  are  quite  elegant  in  their  approach  both  physically  and 
mathematically,  their  predictions  cannot  be  compared  directly  with  experiments 
since  no  relationship  between  experimental  ignition  criterion  and  mathematical 
criteria  has  been  found.  To  make  the  present  gas  phase  distribution  is  not 
known. 

If  at  any  instant,  we  consider  an  elementary  volume  *dv'  in  the  gas 
phase  at  a  distance  x  from  the  propellant  surface,  the  intensity  of  electro- 
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magnetic  radiation  emitted  by  the  species  emitting  radiation  at  a  frequency  v 
is  given  by9; 


die  NrnthV) 

2.18 


,  xi  xt  -Eu/kT 

where  N  =  N  (— )e 
,  m  mo  Q 
and  ^ 


is  the  number  of  molecules  in  the  excited  state  and 


A  = 


mn 


64t t3 


3 

w 


mn 


3h 


is  the  Einstein  Coefficient  of  spontaneous  emission.  Absorption  is  neglected 
since  it  is  generally  compensated  by  the  induced  emission.  The  total  number  of 
molecules  of  species  emitting  radiation  is  related  to  the  number  of  product 
molecular  present  as  NmQ  =  N  f  =  11(1  -  y^  -  y^  -  y  )f  where  N  is  the  total 
number  of  molecular  present  £er  unit  volume  of  the  °XgaI  at  temperature  T. 
Noting  that  dv  =  a  dx  and  defining  L*  as  the  characteristic  slit  width  the 
expression  for  the  intensity  of  light  passing  through  the  width  becomes 

1=  K  Vi)ex|»(-h\>/kT) doc  2  19 

© 


.  „  64  4 

where  K  =  —  tt 

the  shock  tube. 


2 


ac  being  the  cross  sectional  area  of 


This  energy  impinges  on  a  photocell  through  an  aperture  with  area 
'a  '  at  a  distance  ’  r’  from  the  centre  of  the  shock  tube  through  a  suitable 

energy 
-e2/kT 

vp™  k  2.20 

o 

X  0  (x  t )  *—“ 

where  <p  =  /  p? —  dx  and  K  =  N^K.  N^  is  the  number  of  gas  molecules  present 

in  the  initial  phase  per  unitvolume.  The  relative  response  of  the  photo¬ 
detector  is 

J(t.)  *  2.2i 

On  the  basis  of  light  emission,  the  ignition  delay  is  defined  to 
be  the  time  t*,  required  for  J  to  reach  some  value  J*.  phus  the  criterion 
is  the  achievement  of  J  =  J*. 


filter  passing  a  very  narrow  region  of  the  specturm.  Assuming  that  the 

of  all  the  excited  states  are  sufficiently  high,  Q  =  g  +  g  e'Ei/W  , 

.  _£  /kT  °  1  ~ 

+  8je  J  +  ...  =  gQ,  the  voltage  output  of  the  photocell  is  given  Dy 


As  a  check  on  the  applicability  of  this  criterion  in  terms  of 
capabilities  of  commonly  used  photo-detectors,  the  theoretical  intensity  of 
thermal  light  emission  in  the  UV  region  was  estimated.  It  was  found  to  be 
of  detectable  magnitude.  This  was  done  for  a  500p  thickness  region  of  gas 
adjacent  to  the  solid  at  25  atms.  and  3000  K,  containing  approximately  10^Upro- 
duct  molecules.  We  assume  that  the  radiation  is  due  to  characteristic  strong 
CN  emission  bands  in  the  UV  region  involving  (Reference  10)  transitions  at 
3883  X  and  3590  %.  The  intensity  of  radiation  falling  on  a  detector  given 
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by  the  right  hand  side  of  equation  2.19  assuming  —  is  on  the  order  of  10  to 

10"  ^  lumens.  For  an  RCA  type  935  photocell1  ^  preceded  by  a  KODAK  Wratten 
filter  18A  ^  intensities  on  the  order  of  10  lumens  are  easily  detectable. 
Therefore,  it  is  concluded  that  thermally  excited  emission  must  comprise  a 
part  of  the  light  detected  in  shock  tube  experiments. 

It  is  possible  that  a  part  of  the  light  detected  by  the  photo-cell 
may  be  of  chemiluminiscent  nature.  Since  chemiluminiscence  is  attributed  to 
radical-radical  interactions,  its  dependence  on  gas  phase  properties  is  not 
known.  One  could  argue  that  such  interactions  will  not  be  frequent  until 
relatively  large  concentrations  exist,  such  as  after  the  establishment  of  an 
intense,  thick  reaction  zone  involving  external  gas  phase  oxidizer.  Such 
conditions  would  correspond  to  the  'long  time"  behaviour  of  the  gas  phase 
ignition  process.  Thus  it  can  be  argued  that  the  thermally  excited  emission 
used  in  the  present  investigations  is  an  adequate  representation  of  reality. 

The  assumption  was  that  CN  emission  is  of  thermal  origin.  There 
is  some  experimental  evidence  that  this  is  so,  as  is  discussed  in  reference  20. 
However,  till  conclusive  experiments  are  performed,  it  is  not  possible  to  say 
with  any  certainty  that  CN  emission  is  thermal  or  not.  Till  that  time,  it  is 
assumed  that  the  present  analysis  is  valid.  In  case  it  becomes  clear  that 
CN  emissions  are  not  thermal,  it  is  possible  to  consider  the  OH  radical,  which 
is  known  to  undergo  thermal  excitations.  The  results  obtained  using  this 
radical  will  not  be  very  different  from  those  of  CN,  since  the  characteristic 
parameters  defining  the  emission  have  similar  magnitudes. 

The  long  time  criterion  discussed  previously  and  the  strong  emission 
criterion8are  somewhat  related.  Detailed  analysis  of  this  is  given  in 
reference0.  This  gives  some  basis  for  the  assertion  previously  made  that 
detectable  or  strong  emission  occurs  close  to  steady  state. 
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CHAPTER  III 


NUMERICAL  METHODS: 


The  set  of  differential  equations  2.6  to  2.7  with  their  initial  and 
boundary  conditions  are  highly  non-linear  in  nature  and  cannot  be  solved  exactly. 
This  leaves  the  on.  /  other  alternative  of  solutions  by  difference  approximations. 
There  are  several  different  schemes  available  and  the  most  common  ones  are  (1) 
explicit  method  (2)  implicit  method  and  (3)  implicitly-explicit!3  method.  There 
are  other  methods  also  in  use,  like  predictor-corrector  methods,  Crank-Nicolson 
scheme,  collocation  method,  etc.  These  methods  have  been  found  to  be  of  not 
much  value  in  the  present  case  because  of  the  large  number  of  equations  to  be 
solved  and  their  highly  non-linear  nature.  The  discussion  of  several  numerical 
methods  was  done  in  reference  (5)  and  only  a  discussion  of  implicltly-explicit 
method  is  presented  here. 


3.2  The  Implicitly-Explicit  Method: 

In  this  method,  one  set  of  nodes  is  computed  using  a  simple  explicit 
scheme  and  alternate  nodes  are  computed  using  an  explicit  scheme.  For  example, 
consider  a  differential  equation  of  the  type: 


aflUx  +  3.1 

In  this  case,  all  the  eve.,  nodes  can  be  calculated  using  the  explicit 
scheme.  Thus 

K+l  ,  . K  *  k  .  fc  k  ,»  K  k. 

___  ,  - - - +  a<>  ——  4 r21 


where  (i)  refers  to  the  space  point  and  (k)  to  the  time  step.  At  and  Ax  are 
increments  in  time  and  space  variables  respectively.  To  compute  intermediate 
odd  nodes  we  use  the  implicit  difference  approximation 


a 


w  /i*  o.i  **  fiK*'  u*+t  *+> 

+  (X  ~  + 

At  AX.1  ° 


1AX. 


3.3 


The  value  of  f 


k+1 

2i-l 


can  be  approximated  by 


M ) 


3.4 


Since  the  equation  3.2  is  explicit  equation  3.4  is  also  explicit. 
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It  can  be  seen  from  above  that  an  implicitly-explicit  method  is 
basically  an  explicit  method  involving  the  same  number  of  arithmetic 
operations  as  a  simple  explicit  method.  The  advantage  of  an  implicitly- 
explicit  methos  is  obvious.  Since  an  implicit  method  introduces  positive 
error,  i.e.  consistent  overestimation,  and  an  explicit  method  constantly 
underestimates,  the  net  error  is  greatly  minimized.  Hence,  this  method  is 
more  accurate  for  the  same  lumber  of  arithmetic  operations  involved. 
Alternatively,  for  the  same  accuracy  as  either  implicit  or  explicit  method 
alone,  thgfcnumber  of  operations  are  greatly  reduced.  Furthermore,  the  value 

of  A= - 2  which  defines  stability  of  the  approximation  is  unity  instead  of 

(Ax)  ^  as  in  the  case  of  an  explicit  scheme.  Since  the  stability 
criterion  X  =  1  was  deduced  originally-^  for  linear  parabolic  differential 
equations,  for  non-linear  cases  caution  should  be  excercised  to  see  that  X  <  1. 
Convergence  of  the  system  has  to  be  determined  by  trial  and  error  procedure 
as  are  generally  done  for  non-linear  differential  equations. 
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CHAPTER  IV 


DISCUSSION  OF  RESULTS: 

Results  on  homogeneous  propellant  ignition  are  presented  under  three 
headings:  (1)  short  time  results,  (2)  long  time  results  and  (3)  emission  results. 
Short  time  results  correspond  to  weak  criterion  on  temperature,  which  is  the 
establishment  of  a  gas  phase  temperature  50%  in  excess  of  the  initial  gas  phase 
temperature.  Long  time  results  are  related  to  the  attainment  of  quasi-steady 
gas  phase  conditions  and  emission  results  correspond  to  the  achievement  of  a 
certain  specified  intensity  of  thermal  radiation  emitted  by  igniting  gases 
gases. 


4.1  Short  Time  Results : 


At  short  times,  the  ignition  behaviour  can  be  deduced  analytically  by 
some  simplifying  assumptions.  The  validity  of  some  of  these  assumptions  can  be 
justified  from  excellent  agreement  in  behaviour  between  the  numerical  solutions 
and  asymptotic  solutions. 

Figure  (3)  shows  a  logarithmic  plot  of  dimensionless  ignition  delay 
versus  the  Damkohler  number  b,  for  both  feed  back  and  non-  feed  back  cases. 

For  large  initial  pyrolysis  rates,  i.e.  low  b  values,  both  feed  back  and  non¬ 
feed  back  models  give  almost  identical  results  indicating  that  feed  back  is 
not  important  under  these  conditions.  At  low  pyrolysis  rates,  however,  feed¬ 
back  becomes  very  important.  Without  feed-back,  the  gas  phase  will  be 
starved  of  fuel  supply  and  ignition  delay  tends  to  become  extremely  large  as 
indicated  by  sharp  upturns  of  dashed  curves  of  Figure  3.  In  the  feedback 
model,  the  surface  temperature  keeps  on  rising  till  the  fuel  pyrolysis  rate 
is  sufficient  to  cause  rapid  ignition.  At  high  initial  pyrolysis  rates  both 
feedback  and  non  feed  back  models  give  -2/3  slope  in(Figure  3)which  is  in 
agreement  with  asymptotic  prediction.  At  large  y»  values,  the  stoichiometric 
ratio  (n)  has  virtually  no  effect  on  the  ignition  behaviour.  Only  at  low  y~x 
when  the  evolved  oxidizer  plays  an  important  role  in  causing  ignition  does 
the  stoichiometric  ratio  (n)  become  effective.  This  can  be  seen  by  comparing 
Figures  3  and  4. 

It  is  also  evident  from  Figures  (4)  and  (5)  that  pure  fuels  ignite 
faster  than  propellants  which  ir  in  variance  with  experimental  results.  Also 
in  contradiction  with  experiments  is  the  theoretically  predicted  weak  dependence 
of  ignition  delay  upon  y®  .  Employing  short  time  criteria,  both  numerical  and 
asymptotic  results  show  0  that  t* 'v  (ygx)“^/3  ag  shown  in  Figure  10,  curves  3  and  4. 
Experiments  on  the  other  hand  indicate  that  (y®^) 4  where  6  ranges  from  -1 
to  -3.  The  findings  of  the  present  investigations  employing  short  time  criterion 
show  the  same  basic  trends  as  in  reference  (4)  and  (5). 

4.2  Long  Time  Results: 

Continued  numerical  integrations  —  beyond  the  achievement  of  weak 
ignition  criterion  —  of  system  of  partial  differential  equations  presented 
in  Chapter  II  reveals  very  interesting  fact  that  quasi-steady  burning 
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conditions  prevail  soon  after  the  weak  temperature  criterion  i3  satisfied. 

In  the  period  between  the  attainment  of  quasi-steady  conditions  and  the 
achievement  of  weak  temperature  criterion,  the  conditions  in  the  gas  phase 
change  continuously  and  rapidly.  Figure  (6)  shows  gas  phase  distributions 
at  the  time  when  weak  temperature  criterion  was  satisfied.  Oxidizer  mole 
fractions  are  nearly  constant  everywhere  and  equals  to  the  free  stream 
oxidizer  mole  fractions.  This  strengthens  the  original  hypothesis  that  at 
short  times  reactant  consumption  can  be  neglected.  Also  can  be  seen  from 
the  Figure  6  is  a  thin  intense  reaction  zone  adjacent  to  solid  surface.  This 
later  on  shall  be  termed,  for  convenience,  as  primary  reaction  zone.  Primary 
zone  is  a  pre-mixed  reaction  zone.  Figures  (7)  and  (8)  are  schematic 
representations  of  gas  phase  distributions  soon  after  and  a  long  time  after 
the  short  time  temperature  criterion  is  satisfied.  Between  t  =  0.5  milliseconds 
and  t  -  2.0  milliseconds,  conditions  in  the  gas  phase  are  changing  slowly. 
Conditions  in  the  primary  zone  are  nearly  steady.  The  secondary  zone,  so- 
called  because  of  its  structure  is  a  diffusion  reaction  zone  similar  to  the 
diffusion  reaction  zone  in  the  flames  studied  by  Burke- Schumann.  The  secondary 
reaction  zone  is  thick  compared  to  the  primary  zone.  The  fuel  which  is  not 
consumed  in  the  primary  zone  is  used  up  in  the  secondary  zone.  Another  noticeable 
feature  of  the  secondary  zone  is  that  it  is  steadily  moving  away  from  the 
propellant  surface.  It  has  been  found  from  numerical  computations  that  this 
type  of  behaviour  is  established  for  a  wide  range  of  initial  pyrolysis  rates 
of  the  solid  composition  of  the  initial  gas  phase  and  a  range  of  stoichiometric 
(2<n<6)  of  the  gas  phase  chemical  reactions. 

Situations  similar  to  that  revealed  in  Figures  (7)  and  (8)  are  more 
appropriately  described  by  the  term  'quasi-steady*.  The  quasi-steady  state 
requiring  an  order  of  magnitude  greater  in  time  than  the  weak  ignition  criterion 
that  0*  =  1.50  ,  is  still  achieved  rapidly  with  ignition  times  on  the  order  of 
the  few  milliseconds,  even  in  a  neutral  gas  phase.  The  behaviour  of  the  model 
in  this  time  domain  is  termed  as  ’-ong  time'  behaviour.  An  interesting  result 
occurring  in  the  long  time  domain  is  the  formation  of  twin  flame  zone  structure 
in  the  gas  phase  whenever  conditions  permit.  The  secondary  zone,  as  explained 
earlier,  is  formed  between  the  hot  fuel  and  the  external  oxidizer,  the 
strength  of  which  depends  upon:  (1)  the  mole  fraction  of  the  oxidizer  present 
in  the  initial  gas  phase,  (2)  the  stoichiometric  ratio  (3)  propellant  composition. 
A  strong  diffusion  reaction  zone  exists  for  sufficiently  large  values  of  a  and 
n.  The  strength  of  the  reaction  zone  is  measured  by  the  amount  of  heat  released 
in  the  reaction  zone.  At  a  =  1,  i.e.  for  pure  fuel,  only  secondary  zone  exists 
and  for  y~x  =  0  only  primary  zone  exists. 

As  in  all  previous  models,  the  apparent  ignition  characteristics 
exhibited  by  the  present  model  in  the  long  time  domain  depend  upon  the  choice 
of  ignition  criteria.  In  this  domain  several  ignition  criteria  exist,  for 
example;  (1)  time  to  achieve  maximum  absolute  reaction  rate  (2)  time  required 
for  a  zero  oxidizer  mole  fraction  to  occur  at  any  location  (3)  a  specified  low 
rate  of  change  of  heat  released  in  the  secondary  zone,  etc.  The  results  of 
these  choices  are  indicated  in  Figure  (9).  In  determinin  the  closeness  of 
the  system  to  steady  state,  the  important  quantity  is  (—  -— !•) —  i.e.,  the  per¬ 
centage  change  in  heat  release  rate  in  the  secondary  zoSe.This  quantity  may  be 
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represented  by 
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Curves  5  and  6  of  Figure  (9)  illustrate  the  ignition  behaviours  when  a  specified 
value  of 
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was  chosen.  The  most  satisfactory  criterion  was  — - — —  =  e,  where  e  can 


have  any  specified  value.  This  criterion  tends  to  compensate  for  small  values  of 
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when  y»  is  on  the  order  of  0.1.  As  discussed  elsewhere,  this  criterion  is  also 
appropriate. 


Figure  (10)  shows  a  plot  of  t*  versus  yOT  for  several  initial  pyrolysis 
rates  and  propellant  composition.  As  can  be  seen°Irom  the  graph,  it  is  interesting 
to  note  that  there  is  a  strong  dependence  of  t*  upon  yTO  as  was  observed  during 
shock  experiments.  The  results  are  also  in  agreement  with  the  asymptotic  analysis, 
at  long  times.  It  is  also  demonstrated  from  numerical  results  that  the  ignition 
delay  is  independent  of  initial  pyrolysis  rates,  in  agreement  with  experimental 
results?.  Another  interesting  observation  that  can  be  made  from  Figure  (10)  is 
that  a  pure  fuel  ignites  slower  than  a  propellant  under  similar  conditions,  as 
has  been  found  experimentally. 


4.3  Emission  Results 

The  results  of  present  investigation  emphasize,  still  further,  the  con¬ 
clusions  of  the  previous  investigations  on  the  characteristics  of  gas  phase 
ignition  process  of  solid  materials  -  that  choice  of  ignition  criterion  is  decisive 
in  determining  the  apparent  ignition  behaviour.  The  results  also  demonstrate  the 
importance  of  choice  of  ignition  criteria  in  the  interpretation  and  reporting  of 
experimental  data.  r1  ■ 

Figures  12(a)  and  12(b)  are  plots  of  calculated  relative  response  of  a 
hypothetical  photo-detector,  J,  versus  time  for  two  values  of  'activation  energies 
assuming  equal  reaction  rates  for  both  energy  values  at  a  temperature  of  1800  K. 
Values  of  J  of  0.1  indicate  an  absolute  intensity  on  the  order  of  10~^  lumens  as 
detected  by  the  photo-cell.  These  results  can  either  be  used  to  produce  theoretical 
oscilloscope  traces,  or  achievment  of  a  particular  value  of  J  =  J*  to  obtain  t* 
versus  y~x  plots  for  a  particular  J*  value. 

Figure  13(a)  is  a  linear  plot  of  relative  detector  response  versus  time  for 
y*^  =  1.0,  and  0.5  and  illustrates  the  theoretical  prediction  of  the  way  similar 
experimental  curves  should  look  for  a  solid  undergoing  gas  phase  ignition  process. 
Figure  13(b)  is  a  replot  of  experimental  photo-cell  response  versus  time  during  the 
ignition  of  a  composite  solid  propellant  for  y°°  =  1,  and  0,5  exposed  to  a  reflected, 
shock  wave.  As  was  done  in  experiments,  if  we°choose  the  first  discernable  upward 
deflection  or  ’breakup’  for  determining  the  ignition  event,  it  is  obvious  that 
sensitivity  of  the  detecting  equipment  is  important  in  determining  the  instant  of 
ignition.  For  example,  a  decrease  in  sensitivity  by  a  factor  of  3  in  Figure  13(a) 
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(i.e.  a  replot  with  a  vertical  scale  contracted  by  a  factor  of  three)  would 
produce  curves  that  look  very  much  like  those  in  Figure  13(b)  indicating  that 
the  ignition  behaviour  predicted  by  theory  and  observed  from  experiments  may 
change  considerably  depending  upon  the  sensitivity  of  the  detecting  device. 

Figure  14(a)  and  14(b)  are  plots  of  t*  versus  y»  with  ignition 
criteria  of  J*  *»  0.01,  O.O*'  and  0.1.  The  figures  also  include  the  so-called 
break-up  criterion  which  is  the  ignition  delay  obtained  from  theoretical 
oscilloscope  traces.  As  can  be  seen,  the  dependence  of  t*  upon  y»  becomes 
stronger  and  stronger  as  the  value  of  J*  is  increased  from  0.01  t8x0.1. 

Beyond  0.1,  any  further  increase  in  J*  does  not  change  the  dependence  of  t 
upon  y«  very  much,  though  the  magnitude  of  ignition  delay  changes.  Value!8of 
J*  0?1  correspond  to  long  time  criterion  and  J  «*  0.01  correspond  to  weak 
temperature  criterion.  It  can  also  be  seen  from  Figures  14  (a)  and  14  (b)  that 
ignition  data  for  an  E/R  of  6000°K  do  not  differ  much  from  those  for  an  E/R  of 
12000°K.  This  is  expected  since  the  chemical  time  for  both  activation  energies 
was  taken  to  be  the  same  at  1800°K. 

Figure  15(a)  and  15 (b^  allow  comparison  of  the  present  theoretical 
results  using  several  ignition  criteria  and  experimental  results.  For 
description  of  curves  in  Figures  15(a)  and  15(b)  reference  should  be  made  to 
Tables  2(a)  and  2(b)  respectively.  The  long  time  ignition  Curve  5  and  the 
emission  criterion  Curve  4  exhibit  the  same  slope  for  high  initial  oxidizer 
mole  fractions.  Curve  (1)  is  the  data  obtained  at  ONERA  and  Curves  2  and  3 
were  taken  from  reference  (7).  The  slopes  of  Curves  1  and  3  are  in  good 
agreement  with  those  of  4  and  5.  Magnitudes  of  ignition  delays  are  of  the 
same  order  of  magnitude  for  both  theoretical  and  experimental  curves,  except 
for  Curve  (5)  which  is  an  order  of  magnitude  higher  and  probably  would 
correspond  to  ONERA  data.  Curve  1. 

Figure  15(b)  illustrates  a  further  comparison  of  the  present  theoretical 
results  with  those  of  experiments.  Curves  1  and  3  are  experimental,  and  curves 
2  and  4  are  theoretical.  Comparison  of  Figures  15(a)  and  (b)  illustrates,  to 
some  extent,  the  diversity  of  the  experimental  data  existing  for  composite 
propellant  ignition  under  shock  tube  conditions.  They  also  indicate  with 
reference  to  Table  II,  the  effect  of  different  experimental  and  theoretical 
ignition  criteria  on  the  apparent  ignition  behaviour. 

Figure  (16)  illustrates  the  difference  between  the  behaviour  of 
ignition  of  pure  fuels  and  a  homogeneous  propellant  when  an  emission  criterion 
is  employed.  It  can  be  seen  that  the  ignition  of  a  pure  fuel  is  far  more 
sensitive  to  the  initial  oxidizer  mole  fraction  in  the  gas  phase  than  is  the 
propellant.  At  large  y00  values  a  pure  fuel  ignites  faster  than  a  propellant 
for  equal  J*  values  and°lhe  situation  reverses  as  y®  is  decreased.  The  long 
time  ignition  data  shows  that  a  pure  fuel  ignites  slower  than  a  propellant  at 
larger  values  of  y“  .  This  perhaps  is  due  to  the  fact  that  quasi-steady 
conditions  are  attained  in  the  gas  phase  sooner  for  a  propellant  than  for  a 
pure  fuel.  However,  it  is  not  clear  why  an  emission  criterion  shows  a  different 
behaviour  than  the  long  time  behaviour. 

4.4  Pressure  Dependence  of  Ignition: 

Under  shock  tube  conditions,  an  increase  in  the  initial  gas  pressure 
has  a  profound  effect  on  the  jump  temperature  at  the  propellant  surface.  Since 
the  propellant  pyrolysis  rate  is  an  exponential  function  of  surface  temperature 
it  can  be  expected  that  gas  phase  total  pressure  has  considerable  effect  on 
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ignition  timesj  In  Chapter  II,  asymptotic  analysis  at  short  times  revealed 
that  t*  (Aub  y«  )  '  ,  where  o>  =  S—  .  in  this  case  m  is  the  pyrolvsis  rat 


In  this  case  m  is  the  pyrolysis  rate; 
the  pyrolysis  rate  at  40CTK.  Since  p  D  'v  p, 


at  the  jump  temperature  and  mo  is  0  the  pyrolysis  rate  at  400 

T  t/’  i(y»V»Vl*to)  >  ^  ^  /' 706,2 )  an<^  A*  $5R/r>Cp&fc 


we  have 


t*f~  "  Cw/n.„'  ►*] 


To  find  the  dependence  of  (m)  upon  the  initial  pressure  we  assume  a  power  law 
relationship  of  the  form:  r? 

222.  -  /  T~  Tzo  ) 
n>0  [  +GO-Ts,J  0 

and  equating  that  at  any  jump  temperature  other  than  400  K,  this  ratio  should  oe 
the  same  as  . 

e^p  (-  epl  P2TS  )  j  expt-  Bp/+ooR  ) 


from  which  it  was  found  that  n  =  2  for  an  E  /R  of  6000°K.  Thus  we  obtain  the 
pressure  dependence  of  ignition  delay  as  t  ^  -v  p“5/3#  Ac  long  times>  however, 
since  t*  ^(An°°2)  it  can  be  expected  that  t*  ^  p"-1-. 


Figure  17(a)  and  (b)  are  theoretical  oscilloscope  traces  for  mass 
fluxes  of  0.05  and  0.005  g/cni  sec. ,  respectively  at  y»  =  1.  The  ignition 
delays  obtained  from  these  graphs  with  J*  =  0.1,  J*  =  i).05  and  the  break-up 
criterion  are  plotted  on  lnt*  versus  lnp  to  give  ignition  delay  as  a  function 
of  pressure,  in  Figures  18(a)  and  (b) .  It  is  obvious  from  these  figures  that 
t*-p  relationship  strongly  depends  upon  the  choice  of  ignition  criteria. 


With  break-up  criterion  (weak  emission)  or  short  time  temperature 
criterion,  numerical  solutions  give  results  in  very  good  agreement  with  short 
time  asymptotic  relationship  deduced  above.  Thus,  from  the  graphs  it  is  found 
that  t*  'v  p--'-*2^  where  as  asymptotic  relationship  gives  t*  p~l*&7.  This 
is  also  in  agreement  with  experimental  observations. 


As  the  value  of  J*  increases  to  0.05  or  0.1,  the  dependence  of  t* 
upon  p,  predicted  by  the  theory,  decreases  approaching  a  limiting  value  of 
t*  'v  for  large  J*  values.  Here  again  large  emission  results  are  in  agree¬ 
ment  with  long  time  asymptotic  analysis.  Thus  it  can  be  said  that  more 
sensitive  device  may  indicate  a  strong  pressure  dependence  whereas  an  insensitive 
device  may  do  otherwise. 


The  response  of  the  predicted  ignition  delay  of  feed  back  model  to 
changes  in  concentration  of  the  oxidizer  in  the  gas  phase  is  shown  in  figure 
18(c).  It  is  important  to  note  that  under  physical  situation  assumed  in  the 
present  analysis,  and  in  end  wall  shock  tube  experiments,  the  real  ignition 
delay  is  not  purely  a  function  of  oxidizer  concentration  alone.  Different 
ignition  delay  oxidizer  concentration  relationships  are  predicted  for  the 
cases  of  (1)  variable  mole  fraction  at  constant  pressure  and  (2)  variable  total 
pressure  with  constant  oxidizer  mole  fraction,  as  shown  in  Figure  18(c). 
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CHAPTER  V 


CONCLUSIONS 


The  present  theory  succeeded  in  explaining  several  experimental 
observations  made  on  both  double  base  and  composite  propellants.  The 
present  model  is  an  improvement  over  the  previous  ones  in  that  it  can 
adequately  represent  a  propellant  igniting  both  in  oxygen  rich  and  neutral 
environments  and  the  theoretical  predictions  can  be  directly  compared  with 
experiments.  This  was  not  possible  in  previous  models.  In  addition,  it 
has  been  possible  to  evaluate  the  validity  of  some  simplifying  assumptions 
that  have  been  made  in  previous  studies  of  this  type. 

One  of  the  most  important  results  of  the  present  investigations  is 
that  completely  transient  conditions,  including  feedback,  do  indeed  allow 
the  solid  itself  to  achieve  quasi-steady  burning  conditions.  It  is  found 
that  feedback  is  operative  over  a  wide  range  of  decomposition  rates  of  the 
solid.  Nevertheless,  feedback  can  be  neglected  for  sufficiently  high 
decomposition  rates  of  the  solid  without  significantly  affecting  the  final 
results. 


For  solids  producing  fuel  and  oxidising  decomposition  products  in 
non-stoichiometric  amounts,  it  has  been  found  that  a  twin  reaction  zone 
structure  exists  in  the  gas  phase  soon  after  the  short  time  criterion  of  the 
form  6*  =  1.56  is  satisfied.  The  size  and  strength  of  the  secondary  zone 
strongly  depends  upon  the  amount  of  oxidizer  initially  present  in  the  gas 
phase.  It  has  been  found  from  numerical  solutions  that  for  a  given  propellant, 
the  size  and  strength  of  the  primary  reaction  zone  is  independent  of  the 
initial  oxidizer  mole  fraction.  This  supports  the  hypothesis  that,  any 
criterion  corresponding  to  the  development  of  the  primary  reaction  zone,  like 
short  time  or  weak  emission  criterion  indeed  should  exhibit  a  weak  dependence 
upon  the  external  oxidizer.  Since  experiments  indicate  a  stronger  dependence 
of  t*  on  y00  it  is  to  be  concluded  that  experimental  detection  does  not 
correspond °?o  the  development  of  the  primary  zone. 

The  reason  for  the  existence  of  secondary  zone  is  discussed  in  the 
previous  chapter.  It  is  only  to  be  mentioned  here  that  since  most  propellants 
are  fuel  rich  and  shock  tube  experiments  are  conducted  at  sufficiently  high 
oxidizer  mole  fractions  one  can  always  expect  a  twin  flame  structure  in  the 
gas  phase.  The  implications  of  this  statement  is  obvious.  It  points  out 
that  any  criterion  related  to  the  development  of  a  strong  secondary  zone 
should  reveal  properties  entirely  different  from  those  revealed  by  weak  criterion. 
Since  the  strength  of  the  secondary  zone  greatly  depends  upon  the  initial 
oxidizer  mole  fraction  it  is  clear  that  ignitionbbehaviour  should  do  likewise. 

The  above  fact  is  demonstrated  very  well  from  the  fact  that  both  long 
time  criterion  and  large  emission  criterion  indicate  the  same  result  that  t* 

'v  (y~x)-2  and  t*  f  (mQ) .  Experiments  show  similar  behaviour. 


24 


The  development  of  emission  criterion  helps  to  compare  theoretical 
results  with  those  of  experiments  under  the  criterion  of  detectable  amount 
of  light  emission  for  the  hot  gases  as  detected  by  a  photo-cell.  The  light 
emission  criterion  was  derived  by  determining  the  time  required  for  a  quantity 
J,  proportional  to  the  output  of  a  hypothetical  photo-cell  detector  to  reach 
some  selected  value  J*.  This  detector  sensed  the  calculated  intensity  of 
thermally  stimulated  light  incident  on  it. 

The  new  criterion  can  be  considered  as  a  reasonable  theoretical 
equivalent  of  the  experimental  criterion  employed  during  shock  tube  experiments 
and  helps  to  clarify  the  diversity  present  in  the  experimental  data  and  type  of 
ignition  data  needed  to  ascertain  the  ignition  mechanism.  The  theory  based  on 
emission  criterion  predicts  that  a  highly  sensitive  device  should  indicate  a 
weak  ignition  dependence  on  external  oxidizer  mole  fraction  and  a  less  sensitive 
device  should  do  so  otherwise.  Hence  the  theory  calls  for  experiments  using 
at  least  two  greatly  different  sensitivities  of  the  photo  detector  to  determine 
the  mechanism  of  ignition. 

It  has  also  been  found  that  the  secondary  zone  is  responsible  for 
the  majority  of  the  light  detected.  The  weak  emission  criterion  gives 
essentially  the  same  results  as  with  short  time  temperature  criterion  and  the 
strong  emission  criterion  corresponds  to  long  time  criterion.  This  is  quite 
well  substantiated  in  Chapter  II. 

It  was  shown  that  purely  thermal  excitation  of  certain  radical  was 
of  sufficient  magnitude  to  produce  light  intensities  detectable  by  an  ordinary 
photo-cell  similar  to  RCA  935.  It  is  possible  that  some  of  the  light  detected 
experimentally  may  be  due  to  chemi-luminiscence.  Present  analysis  ignores  this. 
Thermally  excited  light  emission  used  in  the  present  investigation  can  be 
considered  as  an  adequate  respresentation  of  reality.  Thus  it  can  be  concluded 
that  ignition  of  fuels  and  homogeneous  propellants  under  shock  pulse  cause 
explained  by  pure  gas  phase  mechanisms . 

In  the  above  discussion  it  was  assumed  that  ignition  was  triggered  by 
pure  gas  phase  mechanisms.  However,  in  reference  (6)  it  is  shown  that  both 
simple  heterogeneous,  and  simple  gas  phase  theories  predict  nearly  the  same 
ignition  behaviour.  From  this  point  of  view,  the  theoretical  analysis  carried 
out  so  far  is  not  really  conclusive  as  to  the  exact  mechanism  of  ignition  during 
shock  tube  tests,  since  equivalent  work  has  not  been  carried  out  for  hetero¬ 
geneous  model.  In  reality,  it  may  be  possible  that  both  mechanisms  play  their 
part  in  the  ignition  process.  The  purpose  of  the  present  investigations  is 
therefore  not  to  invalidate  the  surface  reaction  theory  but  to  show  the 
characteristic  behaviour  of  purely  gas  phase  ignition  process,  and  such  a 
mechanism  can,  in  fact,  explain  all  the  experimentally  observed  ignition 
behaviour  —  something  which  has  not  been  demonstrated  previously. 


25 


CHAPTFR  VI 


6.1  Introduction: 


In  all  the  previous  theories  on  gas  phase  ignition  under  shock  tube 
conditions  it  was  assumed  that  either  the  solid  was  a  pure  fuel  or  contained 
an  oxidizer  which  had  the  same  physico-^J  _j.c  properties  as  the  external 
oxygen.  Thus  it  was  found  that  tbr  aependence  of  ignition  times  y°°  became 
weaker  as  the  amount  of  oxidizes  present  in  the  propellant  was  increased. 

If  the  evolved  oxidizer  is  less  active  than  the  ambient  oxidizer,  it  is 
obvious  that  the  ambient  oxidizer  is  the  controlling  factor  in  deciding  the 
ignition  characteristics.  Therefore,  previous  analysis  applicable  to  pure 
fuels,  discussed  in  Chapter  II,  can  be  employed  in  estimating  the  ignition 
behaviour  of  the  propellants.  On  the  other  hand  it  would  be  interesting  to 
evaluate  the  situation  in  which  the  evolved  oxidizer  is  far  more  active 
than  the  ambient  oxidizer.  In  such  a  case  one  would  anticipate  that  ignition 
behaviour  should  less  strongly  depend  upon  the  external  oxidizer.  An  example 
of  a  propellant  which  produces  an  oxidizer  more  reactive  than  the  ambient 
oxidizer  is  an  ammonium  perchlorate  based  propellant.  This  type  of  propellant, 
on  decomposition  produces  an  oxidant,  perchloric  acid,  which  is  very  active. 
Hence  it  is  of  interest  to  find  out  what  happens  in  such  a  case.  This  is 
discussed  in  the  following  sections. 


6.2  Description  of  the  Model: 


The  assumptions  made  here  are  the  same  in  Chapter  II  except  that  the 
two  oxidizer  diffusion  processes  are  treated  separately  and  the  evolved  oxidizer 
is  assumed  to  have  different  kinetic  properties  than  the  ambient  oxidizer. 

With  these  in  mind,  the  equations  for  conservation  of  reactants  and  energy  can  be 
written  as  below: 


Defining  the  operator 


„  -  3  8 

H  =  t—  -bu>  rr  - 
3t  H 


the  equations  after  proper  non-dimensionalization  may  be  written  as 


H(> n  £ (•- 


H  (%'*)  --%%***)■  (-  Ye) 
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HI  liJ'O  6.4 

H  ;»)=  A?4  'C+*QU**i(''  e)J«K-  '/«)  6.5 

and 

<!&  .  <j6s  _ 

TT  +  s*  75,  ~  Tff  6,6 

Subject  to  the  initial  and  boundary  conditions 


^  (<»•?)  ~o 

n,  (o,T )-  CO-'  ^  ft 
<J* 

*&(**) 

r\LM  i^(0>r;s  6-0  b 

7,(m) 

=  ©(o*,Tjr  e0 

0t(o,x)  *  0(o,r) 

*  eiC"0,,x^-e-f<0 

z54  as,  1  3F 

For  simplicity  it  was  assumed  that  A  =  1,  E=Q=n=l,  but  Z.,  ^  Z_.  All  the 
dimensionless  groups  are  defined  as  in  Chapter  II,  except  for  the  additional 
ones  which  are  defined  as  E  =  E^/E2>  h  =  Z^/Z^t  Q  =  Q^/C^  anc* 


n 
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6 . 3  Approximate  Solutions : 

Limiting  solutions  can  be  obtained  either  for  large  or  small  values  of  v. 
Case  I,  Small  tt: 

Equations  6.1  to  6.6  simplify  to 


H  (>).',  )•  "(1,)*° 

H(V)  •  exK-'/e) 

me)  .  AJitr£e*Kr  '/e> 

i  d@s  .  r /ffO  -j 


The  boundary  conditions  remain  unaltered.  ^  procedure  similar  to  pure  fuel 
case  lead  to  the  relationship  that  x*  (Ab  con,y ”x>~2'3  for  short  times  and 
x*  ^  (An00  )“-*-for  long  times. 

OX 

Case  2,  Large  n 

For  large  II,  equations  6,1  to  6.6  reduce  to: 

- ~ns  exfriO- 

» (%;;  =  "  h  (.10 

h  J  *  -  1hr&  «*K-  ’/e) 

Hie)  =  a  i  )•  '7® ) 
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and 


d6s  f-itu}  &S  _  ^  §£ 

it  +  x/l  <**/ 


With  the  initial  and  boundary  conditions  same  as  before,  the  solution 
in  this  case  is  similar  to  the  case  of  a  propellant  igniting  in  neutral  atmosphere; 
thus  the  ignition  delay  should  be  almost  independent  of  external  oxidizer  mole 
fraction. 

6.4  Results  and  Conclusions: 


In  this  section,  only  the  emission  criterion  is  employed.  Further, 
for  the  reasons  mentioned  above,  only  the  case  of  evolved  oxidizer  more 
reactive  than  the  gas  phase  oxidizer  is  considered.  Figure  (18)  is  a  plot  of 
t*  versus  y»  for  Z..  =  5Z2  or  an  evolved  oxidizer  five  times  as  reactive  as  the 
ambient  oxi8¥zer,  at  the  same  initial  gas  phase  temperature.  It  was  found  that 
for  a  selected  ignition  criterion  of  the  type  J^J*,  dependence  of  t*  upon  y» 
decreases  as  the  reactivity  of  the  evolved  oxidizer  increases.  To  illustrate 
this  effect.  Figure  (18)  contains  two  curves  for  a  =  0.5  and  El/R  =  6000°K  for 
the  cases  of  Z  =  Z  and  Z  =  5Z^.  It  is  evident  that  for  J*  =  0.1,  the  slope 
of  the  curve  for  Z^  =  5Z,  is  lower  than  that  for  Z^  =  Z^.  However,  steeper 
slopes  can  be  expected  at  large  J*  values  since  tne  fuel  that  is  not  used  up  in 
the  primary  zone  forms  a  strong  secondary  zone  with  the  external  oxidizer. 

The  snake  like  character  of  the  lowest  curve  in  Figure  (19)  resulted 
from  making  a  theoretical  oscilloscope  trace  similar  to  Figure  13(a)  and  applying 
a  break-up  criterion  for  the  case  of  Z^  =  5  Z^.  This  probably  is  due  to  the 
difficulty  of  ascertaining  the  exact  break-up  point  from  the  theoretical  oscilloscope 
traces.  It  is  interesting  that  it  bears  a  marked  resemblence .to. composite  pro¬ 
pellant  ignition  data  recently  presented  by  Summerfield  et  alv  .  Such  ignition 
characteristics  may  be  interpreted  as  indicating  a  gas  phase  ignition  process  in 
which  the  evolved  oxidizer  is  considerably  more  reactive  than  the  oxidizer 
initially  present  in  the  gas  phase. 

Since  shock  tube  experiments  in  composite  propellants  indicate  a  large 
effect  of  y»  on  t*,  the  gas  phase  theory  currently  formulated  leads  to  the 
conclusion  ?ftat  either  little  or  no  oxidizer  was  decomposed  during  the  primary 
ignition  process  or  the  evolved  oxidizing  species  have  kinetic  properties  similar 
to  oxygen  during  shock  tube  experiments. 
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CHAPTER  VII 


GAS  PHASE  IGNITION  OF  HETEROGENEOUS  SOLID  PROPELLANTS 
UNDER  SHOCK  TUBE  CONDITIONS 


7.1  Introduction: 

The  present  day  propellants  are  heterogeneous  in  nature.  They 
consist  of  oxidizer  crystals,  like  ammonium  or  potassium  perchlorate,  finely 
ground  and  dispersed  in  a  fuel  matrix.  Upon  decomposition  of  such  a  pro¬ 
pellant,  the  fuel  oxidizer  vapours  may  diffuse  both  axially  and  radially, 
probably  affecting  the  ignition  behaviour  considerably.  Since  the  size  of 
the  oxidizer  particle  is  of  the  same  order  of  magnitude  as  the  reaction  zone 
thickness  (=  20  to  200  y)  it  can  be  expected  that  radial  diffusion  may  play 
a  significant  role  in  determining  the  overall  ignition  characteristics.  To 
extend  the  analysis  cf  gas  phase  ignition  characteristics,  making  the  model 
more  realistic  and  applicable  to  heterogeneous  propellants,  a  two  dimensional 
model  is  considered,  and  a  detailed  description  of  the  model  is  presented  in 
the  following  sections. 

7.2  Physical  Model: 


A  schematic  representation  of  the  physical  model  and  the  propellant 
geometry  are  shown  in  Figures  (20)  and  (21) .  The  assumptions  made  during  the 
equations  are  the  same  as  in  Chapter  II,  except  for  a  few  additional  ones. 

They  are: 

1.  The  diffusion  in  both  axial  and  radial 

2.  The  decomposition  rates  of  fuel  and  oxidizer  are 
assumed  to  proceed  according  to  an  Arrhenius  law. 

3.  Negligible  amount  of  propellant  is  decomposed  during 
ignition  and  the  propellant  surface  remains  planar  during 
rapid  ignition  process  (although  there  may  exist  a  radial 
surface  temperature,  and  hence  mass  flux,  distribution). 

4.  Density  is  constant  in  the  gas  phase. 

5.  Oxidizer  and  fuel  have  identical  thermal  properties. 

It  should  be  noted  that  although  the  assumption  (3)  is  not  necessary, 
assuming  a  planar  interface  simplifies  the  problem  considerably.  In  principle, 
it  is  possible  to  set  up  boundary  conditions  which  take  into  consideration  the 
variation  of  mass  flux  with  radius  and  the  surface  geometry  acquired  by  the 
propellant  surface  because  of  this.  But  the  computational  labour  and  time" 
involved  did  not  permit  the  consideration  of  a  more  exact  model.  Nevertheless 
it  was  felt  that  the  effect  on  the  ignition  behaviour  of  the  boundary  movement 
would  be  insignificant. 
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subject  to  the  following  initial  and  boundary  conditions: 

§><?,{»«}•,  r>0,  £•  00  j  %s~  -eo 

7*  7«  *  C  ,  7i  =  7* 
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In  the  above  analysis  it  was  assumed  that  the  evolved  oxidizer  had  the 
same  physical  and  kinetic  properties  as  the  free  stream  oxidizer.  The  dimensionless 
quantities  appearing  in  the  above  equations  are  defined  as  follows. 
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For  a  significant  part  of  the  work  it  was  assumed  that  the  decomposition  rates 
of  both  oxidizer  and  fuel  were  the  same  at  equal  surface  temperacures.  The 
above  equations  were  numerically  solved  on  an  IBM  360/75  digital  computer. 

The  results  of  numerical  computations  are  presented  in  Section  7.6. 


7 . 4  Asymptotic  Solutions: 


Solutions  applicable  under  special  conditions  can  be  obtained  by 
neglecting  reactant  consumption  up  to  ignition.  If  we  further  assume  that 
each  fuel  strand  is  surrounded  by  huge  oxidizer  strand  inverse  geometry  of 
the  previous  section) ,  the  equations  discussed  above  can  be  reduced  to 


.  <l*'l *  + Xi.fr  £!») 

'  TT‘  r  •It'-  Ttl 
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with  n  *  n°°  for  a  large  ygx  values.  Wall  temperature  was  assumed  constant  and 
convec2£on  Selected.  Negligible  amount  of  oxidizer  decomposes  during 
ignition. 


Equations  7.7  and  7.8  are  subjected  to  initial  and  boundary 
conditions : 


O  *  fSfa  ; 
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Also  at  t  =  0,  •-  0  for  all  £  and  p. 

Successive  application  of  Hankel  and  Laplace  transform  gives  the 
i  or  ri^  as 


to 


o 


solution 
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If  we  define  £,  p ,  x 


r- 

t  =  (A*£nbf4AH 

<j>  =  r^: 


then,  substitution  of  7.9  in  7.8  gives 


7.10 


7.11 


where 


ys 
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Solution  of  7.‘-  J  of  the  form  0  = 0 (t ,5 ,p ,p £  ^x  )  and  thus  for  an 

ignition  criterion  of  the  i_>pe  L  '(0)=  0*  occurring  at  x  =  x*  we  get  x*  'v  1/(A 

n”  nbp2)or  t*  My00..)-1.  This  is  an  interesting  result,  indicating  the 

ignition  delay  depends  inversely  on  the  oxidizer  mole  fraction,  which  is  in 

good  agreement  with  some  experiments.  The  present  result  that  t*  ^  is 

an  interesting  contrast  to  the  previous  result  for  homogeneous  propellants  of 

t*  'v,  (y«  ) — 2/3  and  is  a  direct  result  of  the  2-diraensional  nature  of  the 
i-i  ox 
problem. 


The  above  result  is  equally  valid  for  a  fuel  igniting  in  a  shock 
tube,  rich  in  oxygen,  if  the  shock  tube  diameter  is  much  larger  than  the 
propellant  strand,  as  is  the  case  in  the  experiments  reported. 


7.5  Numerical  Methods: 

Both  explicit  and  implicit  schemes  discussed  in  Chapter  III  can  be 
employed  in  the  present  situation  also.  However,  it  is  obvious  that  either 
of  the  method  will  take  extremely  long  computational  times.  For  example,  for 
the  same  number  of  axial  nodes  and  the  same  number  of  time  steps,  a  2-D  problem 
will  take  roughly  M  times  the  computational  time  needed  by  a  1-D  problem,  where 
M  is  the  number  of  nodes  in  the  radial  direction.  The  so-called  alternating 
implicit-explicit  method  would  seem  to  be  promising,  and  a  rough  estimation 
indicates  that  this  method  is  several  times  as  fast  as  either  the  implicit  or 
explicit  method.  However,  it  is  quite  cumbersome  and  hence  only  the  explicit 
method  is  considered  here. 
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Consider  a  differential  equation  of  the  type 

*  Uxx-t  Uy,+  &+  7.12 

At  all  grid  points  except  at  r  =  0,  and  r  =  R  the  equations  can  be  approximated 
by  the  forward  difference  approximation: 


M  K 
U(j -  Ufj 


i 
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where  x  =  i  h. 


=  jh„  and  t  -  X,k. 


The  above  scheme  offers  no  computational  problem  except  at  r  =  0. 
At  r  =  0,  a  singularity  exists  in  the  term  I  |H.  Since  it  was  assumed  that 


at  r  =  0,  ^  =  0  because  of  symmetry. 


term 


is  of  the  form  0/0 


at  r  =  0  which  is  indeterminate.  If  however,  u* (r)  isrexpanded  in  McLaurin 
series  and  the  limiting  process  is  carried  out  by  application  of  L'Hopital's 
rule,  one  obtains 


to  f  *(&]■“’ W 


7.14 


Hence  along  the  line  r  =  0,  the  equation  to  be  solved  is: 

in.  ,  +  t  +  f  7.15 

It  di*  ir‘ 


At  the  outer  boundary  r  =  R 


and 
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0  and  therefore  equation  becomes 


c 
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With  these  precautions  in  mind,  equation  7.12  can  be  solved  in  a  straight 
forward  manner. 


Stability  and  Convergence: 

Though  rigorous  proofs  of  stability  and  convergence  are  not  possible 
for  the  above  equations,  many  times  stability  of  the  approximation  is  assumed. 
For  a  two-dimensional  linear  problem  it  can  be  shown  that  all  explicit  schemes 
are  stable  as  long  as  ^  <.  Vhf  +  l/h\ 

If  we  choose  that  h^  “  1^  =  h,  then  the  restriction  on  the  grid  size  for 
stability  becomes ^  ifhr  1/4,  a  more  stringent  restriction  than  that  for  a 
one  dimensional  case. 


As  far  as  convergence  is  concerned,  it  can  only  be  tested  by  a  trial 
and  error  procedure.  If  the  grid  size  Is  reduced  to  half  of  its  original  size 
and  the  solutions  thus  obtained  do  not  differ  very  much  from  each  other  then 
the  difference  approximation  can  be  considered  convergent.  Computer  trials 
indicated  this  t<<  be  true  in  the  present  instance. 


7.6  Ignition  Criteria: 


Similar  to  homogeneous  propellant  ignition  data,  data  for  heterogeneous 
propellants  are  presented  using  short  time  and  emission  criteria.  The  form  of 
expression  for  the  intensity  of  emitted  radiation  in  terms  of  output  of  photo¬ 
cell  is  modified  to  p 

Vp"  7!  j  jO-VWi)MK'’  7.17 

1  0  0 

where  p^  is  tne  outer  diameter  of  one  compositt  microstrand.  For  presenting 
the  theoretical  ignition  data,  the  quantity  calculated  during  numerical  com¬ 
putation  was 

J1  ~  X  Z.jfWie“(S  (•  gj-)  7.18 

i-J*o  J 


which  is  related  to  the  homogeneous  propellant  ignition  criterion  J  as 


J  « 


a  J 


where  M  is  the  number  of  nodes  in  the  radial  direction,  assuming  that  the 
grid  size  in  the  axial  direction  is  the  same  for  both  homogeneous  and 
heterogeneous  propellants. 
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CHAPTER  VIII 


8.1  Discussion  of  Results; 

Ignition  behaviour  of  a  heterogeneous  propellant  differs  from  that  of 
a  homogeneous  propellant  in  many  ways.  Firstly,  as  can  be  expected,  there  Is 
diffusion  in  both  axial  and  radial  directions  and  thus  the  degree  of  chemical 
reaction  varies  from  point  to  point.  Hence  it  is  not  possible  to  talk  of  one  or 
two  discrete  reaction  zones  as  was  possible  in  the  case  of  a  homogeneous  pro¬ 
pellant.  Secondly,  since  the  surface  temperature  of  the  propellant  varies 
radially,  the  mass  flux  at  the  propellant  surface  is  not  uniform  resulting  in 
widely  differing  rates  of  cross  diffusion  of  fuel  and  oxidizer  vapour  as  a 
function  of  time. 

Figure  (22)  is  a  plot  of  reactant  mole  fraction  versus  dimensionless 
distance, £  ,  at  a  time  X  ®  0.153  which  is  only  a  short  time  after  the  achieve¬ 
ment  of  a  gas  phase  temperature  equal  to  1.50  .  As  can  be  seen  from  the  figure, 
close  to  the  outer  radius  of  the  propellant  strand,  the  propellant  behaves  like 
a  pure  fuel  since  there  is  more  fuel  at  the  outer  radius  than  at  the  inner  radius. 
Towards  the  centre  of  the  propellant  strand  the  concentration  profile  of  the 
reactant  changes.  It  appears  that  for  small  values  of  r/r  ,  where  r  is  the 
outer  radius  of  the  propellant  strand  (i.e.  close  to  the  centre  of  tRe  propellant 
strand) ,  there  is  a  fuel  hill  established  in  between  the  evolved  and  the  external 
oxidizer.  The  fuel  vapours  flow  down  the  hill  on  either  side  to  form  a  double 
reaction  zone  which  reduces  to  a  single  reaction  zone  as  r/r  is  increased.  It 
is  obvious  from  the  figure  that  the  concept  of  a  single  reaction  zone  is  not 
applicable  to  heterogeneous  propellant  igniticr..  Hence  the  long  time  ignition 
criterion  employed  in  the  case  of  homogeneous  propellants  becomes  meaningless. 

Figure  (23)  throws  more  light  on  the  structure  of  the  gas  phase  during 
the  ignition  of  a  heterogeneous  propellant.  The  fuel  concentration  is  maximum 
at  the  outer  radius  and  the  oxidizer  mole  fraction  is  maximum  at  the  centre  all 
the  time.  Further,  the  fuel  mole  fraction  falls  off  very  sharply  with  distance. 
The  distributions  at|£  =  0.753  do  not  differ  very  much  from  those  at  X  ~  0.153 
except  that  as  the  radius  increases,  the  oxidizer  mole  fraction  first  decreases 
and  then  increases. 

Figure  (24)  is  a  plot  of  relative  emission  intensities  versus  time  for 
y°°  of  1,  0.5  and  0.1,  for  a  particle  size  of  70  microns.  In  Figure  25,  Curve  (1) 
is  obtained  from  short  time  data  and  exhibits  a  slope  of  approximately  -0.8  on  the 
log-log  plot  and  is  in  good  agreement  with  the  approximate  solution  discussed  in  t 
the  previous  chapter.  However,  it  should  be  remembered  that  the  approximate 
solutions  are  valid  only  for  large  oxidizer  particles,  i.e.  for  small  values  of 
^2~^1^2*  Curve  (2)  is  obtained  from  emission  characteristics  and  as  can  be  seen 
from  the  figure,  is  in  good  agreement  with  the  experimental  curve  (3)  for  large 
y°°  values.  Curve  (3)  is  extracted  from  reference  (7).  The  slope  of  curve  (4) 
is  not  very  much  different  from  that  of  curve  (2).  Only  the  magnitude  of  the 
ignition  delay  is  higher. 


37 


Figure  (26)  is  again  a  linear  plot  of  J  versus  t,  for  a  particle  size 
of  200  microns,  It  can  6e  seen  from  the  figure  that  even  at  the  ’break-up' 
point,  ignition  delays  are  very  sensitive  to  the  initial  oxidizer  mole  fraction. 
However,  if  we  take  a  certain  specified  deflection  on  the  oscilloscope  to  define 
the  instant  of  ignition,  as  can  be  seen  from  Figure  27,  for  a  J*  value  of  0.05, 
then  the  dependence  of  ignition,  delay  upon  the  initial  oxidizer  mole  fraction 
becomes  weaker.  A  J*  value  of  10""^  on  Figure  (27)  corresponds  to  the  weak 
temperature  criterion  and  exhibits  a  slope  of  almost  -1  on  the  log- log  plot. 

This  justifies  the  validity  of  approximations  made  in  obtaining  the  asymptotic 
solutions.  It  is  obvious  that  consideration  of  inverse  geometry,  i.e.  fuel 
surrounded  by  oxidizer  crystals,  should  essentially  give  the  same  results  as 
the  normal  geometry,  i.e.  oxidizer  particles  surrounded  by  the  fuel  matrix. 
Figure  (27)  again  demonstrates  that  a  highly  sensitive  equipment  detects  the 
ignition  behaviour  to  be  weaker  than  a  low  sensitivity  device. 

Comparison  between  Figures  (25)  and  (27)  shows  that  the  particle 
size  has  a  strong  effect  upon  the  ignition  delay.  For  example,  curve  (2)  of 
Figure  (25),  which  is  obtained  by  applying  the  break-up  criterion  has  a  slope 
of  -1.2  and  the  curve  (2)  of  Figure  (27)  obtained  by  applying  the  same 
criterion  indicates  a  slope  of  -2.0.  It  can  be  expected  from  these  that  as 
the  particle  size  increases,  the  apparent  ignition  behaviour  exhibits  a 
stronger  dependence  upon  the  external  oxidizer.  As  the  particle  size  is 
reduced,  the  propellant  becomes  less  and  less  heterogeneous  in  nature  and  the 
ignition  behaviour  should  essentially  be  the  same  as  predicted  for  homogenous 
propellants. 

From  numerical  computations  it  has  been  found  that  the  emission  from 
a  heterogeneous  propellant  is  lower  than  the  emission  from  a  homogeneous  pro¬ 
pellant.  This  is  to  be  expected  since,  in  the  case  of  a  homogeneous  propellant, 
the  temperature  inside  the  secondary  reaction  zone  is  high  and  uniform  across 
the  section  whereas  in  the  case  of  a  heterogeneous  propellant ,  the  temperature 
is  high  at  fuel-rich  areas  and  low  at  fuel-poor  areas. 

Continued  numerical  integration  of  the  differential  equations 
describing  the  ignition  of  a  heterogeneous  propellant  indicated  that  quasi¬ 
steady  burning  conditions  are  attained  even  for  heterogeneous  propellants. 

It  has  been  found  from  computations  that  emission  intensities  become  constant 
with  time  when  quasi-steady  gas  phase  conditions  are  attained.  Hence  steady 
burning  conditions  can  be  said  to  be  attained  when  the  emission  intensities 
become  constant  with  respect  to  time. 


Pressure  Dependence  of  Ignition  Delay: 

As  discussed  in  Chapter  IV,  Section  4,  a  change  in  the  total  pressure 
results  in  a  corresponding  change  in  the  initial  jump  temperature  at  the  pro¬ 
pellant  surface  and  hence  a  pressure  dependence  similar  to  that  of  homogeneous 
propellants  can  be  expected.  The  asymptotic  solution  of  Section  7.4  indicates 
that  ignition  delay  should  at  least  be  inversely  proportional  to  pressure  at 
high  mass  fluxes  and  large  y”  values.  Figure  28(a)  and  28(b)  are  theoretical 
'oscilloscope'  traces  for  mass  fluxes  of  0.05  and  0.005  g/cm^sec  at  y®  =1  and 
are  used  to  obtain  £nt*  versus  log  p  curves  for  several  ignition  criteria. 

^Figure  (29)  is  a  plot  of  lot  t*  versus  log  p  for  a  mass  flux  of 
0.05  g/cm  -sec  at  y®  =1.0.  As  can  be  seen  from  the  figure,  at  break-up,  the 
dependence  of  ignit2§n  times  upon  the  gas  phase  total  pressure  can  be  read  as 
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*■1  <  5 

t*«y*  P  ,  The  dependence  does  not  seem  to  change  very  much  if  higher  values  of 
J*  are  chosen,  Only  the  magnitude  of  the  ignition  delay  changes. 

At  low  mass  fluxes  the  pressure  dependence  of  ignition  delay  becomes 
much  stronger.  Figure  (30)  is  again  a  plot  of  log  t*  versus  log  p  for  a  mass 
flux  of  0.005  g/cm^sec&y00  3  1.0,  As  can  be  seen  from  the  graph,  the  pressure 
dependence  at  break-up  is  much  stronger  than  at  high  emission  intensities.  The 
dependence  can  be  read  as  t*v*  p~^'^  at  break-up  and  t*v/»  p“2.0  an  ignition 
criterion  of  J*  =  0.1  is  chosen.  These  are  in  excellent  agreement  with  the 
experimental  resultsl.  Since,  after  the  short  time  criterion  is  satisfied,  the 
mass  flux  from  the  propellant  surface  becomes  sufficiently  high  —  irrespective 
of  the  initial  mass  fluxes  due  to  high  surface  temperatures  attained  in  the 
feed  back  model  —  the  initial  mass  flux  should  not  have  much  effect  on  the 
pressure  dependence  if  a  strong  emission  criterion  (J*  >  0.1)  is  employed. 

8.2  Conclusions 


The  present  theoretical  investigations  on  the  ignition  of  heterogeneous 
propellants  have  revealed  some  interesting  aspects.  Firstly,  with  the  present 
model,  it  is  possible  to  explain  the  experimentally  observed  ignition  behaviour 
of  composite  propellants.  Even  with  a  weak  temperature  criterion  one  theoretically 
obtains  slope  of  -1  on  a  log-log  plot  of  t*  versus  y  which  is  in  good  agreement 
with  some  experiments'.  With  more  stringent  ignitionXcriteria  t*i/> (y°* )~^  or 
stronger  is  found  theoretically  as  well  as  in  experiments"*".  The  theory  indicates 
that  a  sensitive  equipment,  detecting  ignition,  should  indicate  a  behaviour  similar 
to  short  time  criteria.  Theoretical  plots  similar  to  experimental  oscilloscope 
traces  reveal  that  a  criterion  similar  to  experimental  criteria  of  ignition  can 
explain  the  strong  sensitivity  of  ignition  behaviour  upon  the  initial  oxidizer 
mole  fraction  and  the  diversity  of  experimental  results,  A  weaker  dependence  of 
ignition  upon  the  external  oxidizer  should  be  observed  for  propellants  with 
finely  ground  oxidizer  particles.  This  is  because  the  heterogeneity  of  the 
propellant  decreases  as  the  particle  size  is  decreased,  and  f  ister  mixing  is 
achieved  in  the  radial  direction. 

In  the  limiting  case  where  the  particle  size  is  extremely  small  the 
system  should  behave  like  a  homogeneous  system.  (This  can  also  be  explained  by 
the  fact  that  the  diffusion  mixing  time  in  the  radial  direction  is  proportional 
to  the  square  of  the  particle  size  and  consequently  a  reduction  in  the 
particle  size  by  50%  reduces  the  diffusional  mixing  time  by  a  factor  of  3  thus 
making  the  system  more  homogeneous).  The  present  model  also  exhibits  a  strong 
dependence  of  ignition  times  upon  external  total  pressure  which  is  in  good 
agreement  with  the  experimental  results  under  shock  tube  conditions.  On  the 
basis  of  the  present  investigations,  it  can  therefore  be  concluded  that  ignition 
of  composite  and  homogeneous  propellants  under  shock  tube  conditions  can  be 
explained  by  pure  gas  phase  mechanisms,  employing  light  emission  criteria. 


8 . 3  Comparison  Between  Homogeneous  and  Heterogeneous  Propellant  Ignition: 

Tt  has  been  found  from  the  present  investigations  that  under  similar 
conditions  both  homogeneous  and  heterogeneous  propellants  behave  similarly.  In 
the  case  of  a  hcmogeneous  propellant  one  can  speak  of  two  distinct  reaction 
zones,  primary  and  secondary  reaction  zones,  whereas  in  the  case  of  the  hetero¬ 
geneous  propellants,  theoretically  speaking,  the  concept  of  distinct  reaction 
zones  is  not  valid.  For  similar  kinetic  and  physical  properties,  emission 
from  a  homogeneous  propellant  is  more  intense  than  from  a  heterogeneous  pro- 
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pellant.  This  is  because  the  temperature  at  any  section  of  the  gas  phase  is 
uniform  for  a  homogeneous  propellant  and  the  temperature  varies  radially  for 
a  heterogeneous  propellant  with  low  temperatures  near  the  oxidizer  rich  regions. 
For  some  sensitive  photo  detector  read-outs,  heterogeneous  propellants  exhibit 
stronger  sensitivity  for  external  oxidizer  than  do  homogeneous  propellants, 
providing  the  oxidizer  particle  size  is  large.  This  situation  disappears  as 
the  particle  size  decreases. 


Ignition  Behaviour  of  a  Low  Oxidant  Propellant: 

It  is  obvious  from  what  has  been  said  so  far  that  a  propellant  contain¬ 
ing  more  than  50%  oxidizer  by  weight  attains  steady  burning  conditions  after  the 
ignition  transients  are  over.  A  propellant  with  less  than  20%  oxidizer  behaves 
like  a  pure  fuel.  Therefore  it  can  be  expected  that  a  propellant  containing 
oxidizer  between  20  and  50%,  should  behave  like  a  pure  fuel  at  low  initial 
oxidizer  mole  fractions  and  like  a  propellant  at  high  initial  oxidizer  mole 
fractions.  This  effect  is  observable  if  the  decomposition  rates  of  oxidizer 
and  fuel  are  assumed  to  be  different. 

Figure  ('SI)  shows  a  lot  t*  plotted  against  log  y*>  for  a  propellant 
containing  50%  oxidizer  and  35%  oxidizer  respectively.  The  particle  size 
(i.e.  diameter  of  oxidizer  cyrstals)  is  50  to  60%  microns.  A  two  dimensional 
model  was  used  to  compute  the  ignition  data.  It  was  assumed  in  the  numerical 
computations  that  the  activation  energy  of  pyrolysis  of  the  oxidizer  particles 
is  twice  that  of  fuex  and  both  fuel  and  oxidizer  have  equal  decomposition  a 
rates  at  400  K.  As  can  be  seen  from  the  figure;  increasing  the  content  of 
oxidizer  in  the  propellant  only  slightly  increases  the  ignition  delay  at  high 
initial  oxidizer  mole  fractions.  For  an  a  of  0.5  as  y°°  is  decreased,  the 
t*  -  y<»  curves  bends  over  and  at  very  low  y°°  values,  ignition  delay  is 
initiaKy  independent  of  y°°  .  For  a  <  0.3  i?Xhas  been  found  that  ignition 
delay  steadily  increases  as  y°°  is  decreased.  At  a  <  0.35,  there  is  a  bend 
over  or  independency  of  t*  on  y<»  over  a  short  range  of  y°°  values  and  for 
sufficiently  low  values  of  y°°  .  The  bend  over  is  expecte8Xto  occur  at  higher 
values  of  yOT  for  a  stronger  emission  criterion,  J*  >  0.005.  However, 
computational  times  are  already  very  large  and  to  attain  J*  values  of  0.05  or  q,i 
may  take  computational  times  on  the  order  of  an  hour  or  so. 

A  behaviour  similar  to  the  above  has  been  found  during  shock  tube 
ignition  of  composite  propellant^),  Thus  it  can  be  said  that  a  two- 
dimensional  gas  phase  model  described  in  the  present  situation  can  explain 
a  variety  of  shock  tube  ignition  data  presented  by  McAlevy  &  Summerfieid''  1' 
and  Kashiwagi^^  . 
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CHAPTER  IX 

RADIANT  IGNITION  OF  HOMOGENEOUS  SOLID  PROPELLANTS 


9.1  Introduction: 

Arc  image  or  radiant  ignition  of  solid  propellants  has  become  a  wide 
spread  practice  for  screening  propellants  according  to  their  ignitability . 

However,  the  validity  of  application  of  such  data  to  the  ignition  of  propellants 
in  practical  rockets  is  questionable.  In  fact,  there  is  as  much  justification 
for  the  application  of  shock  tube  experimental  data,  to  the  operational  motor 
systems  as  there  is  for  the  application  of  arc  image  experimental  data.  Even 
though  direct  application  of  arc  image  or  shock  tube  ignition  data  for  pro¬ 
pellants  igniting  in  practical  rocket  motors  is  not  possible,  still,  some 
understanding  of  the  ignition  behaviour  of  propellants  is  achieved  by  these 
experiments. 

The  importance  or  significance  of  experimental  data  often  times  cannot 
be  appreciated  or  properly  utilized  without  an  adequate  theoretical  model.  Based 
on  certain  experimental  observations,  generally,  a  simple  theoretical  model  is 
proposed  and  the  predictions  of  the  theoretical  model  are  compared  with  experi¬ 
mental  observations.  Any  discrepancy  is  then  corrected  by  refining  the  original 
model  or  by  proposing  a  new  model.  This  has  been  the  situation  in  the  theoretical 
studies  on  ignition  of  solid  propellants. 

As  already  mentioned  in  Chapter  I,  theoretical  studies  on  the  ignition 
of  solid  propellants  subjected  to  radiant  heating  were  initiated  by  Hicks.  His 
original  theoretical  model  remained  untouched  for  quite  some  time,  till  the 
development  of  composite  propellants.  Somehow,  the  predictions  of  Hicks  model 
did  not  agree  with  the  ignition  behaviour  of  composite  propellants.  Many 
experiments  on  composite  propellants  indicated  that  the  original  hypothesis  of 
Hicks  —  that- ignition  was  caused  by  solid  ■'phase  reactions  —  was  not  correct. 

This  led  to  the  emergence  of  gas  phase  and  heterogeneous  ignition  models. 

Definitive  arguments  and  proof,  for  or  against  either  of  the  new  theories,  is 
beyond  the  scope  of  present  work.  It  is  assumed  here  that  ignition  in  arc 
image  experiment  is  due  to  gas  phase  processes;  the  purpose  is  to  provide 
information  on  just  what  kind  of  ignition  behaviour  might  be  expected  if,  indeed, 
the  ignition  process  in  a  real  propellant  were  entirely  due  to  gas  phase  reactions. 

The  gas  phase  ignition  model  developed  by  McAlevy  and  Summer field  to 
describe  the  ignition  behaviour  under  shock  tube  conditions  was  applied 
successfully  by  Shannon  and  DeverallU?)  to  study  the  ignition  behaviour  of 
solids  igniting  in  an  arc  image  furnace.  A  similar  model,  with  some  refine¬ 
ment,  was  later  on  developed  by  Ohlemiller  and  Summerf ield (18) .  Both  these 
models  showed  considerable  qualitative  agreement  with  experimental  observations. 
Discussion  of  Ohlemiller 's  model  is  only  qualitative  and  exact  solutions, 
either  analytical  or  numerical,  have  not  been  presented  as  yet.  The  object  of 
the  present  investigation  is  to  refine  the  model  of  Ohlemiller  to  take  into 
consideration  the  effect  of  absorption  of  radiant  energy  by  gas  phase  reactants 
and  products  and  make  the  model  describe  the  ignition  of  propellants,  in 
general,  rather  than  pure  fuels  only. 

No  orginality  is  involved  in  the  derivation  of  the  equations  and 
therefore  only  their  final  form  is  presented  here.  Employing  the  assumptions 
listed  in  Chapter  II,  the  equations  of  conservation  of  species  are  written 


41 


as  follows: 
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The  in-depth  pyrolys&s  is  given  by; 
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In  the  above  treatment,  photomeric  dissociation  of  the  propellant  is  not  con¬ 
sidered.  All  the  dimensionless  groups  are  defined  as  in  Chapter  II.  1'^  and 
and  A  are  constants  with 


and 


A  =  '’Wm. 


respectively.  Q  is  the  flux  incident  on  the  propellant  surface  and  is  related 
to  the  source  f’ux  as 

where  xm  is  the  mean  distance  in  the  gas  phase  containing  fuel  species  through 
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which  the  radiation  passes  before  impinging  on  the  propellant  surface.  It 
is  approximated  as 


5Cm  ~ 


J  Sfd*/9(e 


where  y  is  the  wall  fuel  mole  fraction.  Radiation  passes  entirely  through 
oxygen^and  nitrogen.  At  later  stages,  when  significant  amount  of  products 
are  formed,  the  absorption  by  products  should  also  be  considered.  However, 
during  radiant  ignition  negligible  amount  of  products  are  formed  till  the 
fuel  is  ignited.  Also,  it  is  generally  assumed  that  remains  constant 

with  respect  to  pressure.  At  constant  pressure  this  is  the  grey  body 
assumption. 


Since,  during  arc  image  experiments,  the  only  source  of  heat  supply 
to  the  gas  phase  is  the  propellant  surface  upon  which  the  radiant  energy  is 
incident,  the  assumption  that  exothermic  chemical  reactions  are  nearly  absent 
till  ignition  occurs  appears  to  be  quite  valid.  However,  at  extremely  high 
rates  of  energy  supply  to  the  grey  propellant  surface,  the  surface  temperature 
attains  a  sufficiently  high  value  very  rapidly;  hence  the  chemical-diffusion 
time,  t  ,  comprises  the  major  portion  of  the  delay  period.  One  can  rest  assured 
that  such  a  situation  will  not  exist  until  the  flux  level  reaches  several 
hundred  calories/cin>sec. ,  when  the  propellant  Is  not  truly  opaque  and  energy 
is  absorbed  in  depth.  In  reference  it  is  shown  that: 
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where  Q  (1-r)  is  that  portion  of  the  incident  flux  which  is  absorbed  by  the 
propellant  andAT*  is  the  characteristic  temperature  rise  at  the  propellant 
surface  at  the  instant  of  ignition.  In  the  general  case,  the  ignition  time 
is  a  combination  of  these  and  hence 
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which  has  shown  by  Ohlemiller  to  be  in  agreement  with  experimental2results. 
Present  numerical  calculations  show  that  even  for  Qq  of  100  cal/cm  sec., 
t  ,<.<  t .  based  on  e*  -  1.5ft,  or  emission  orae*  =  constant,  criterion, 
unless  a^ery  low  value  of  A©*°is  assumed.  A  detailed  discussion  of  the 
numerical  results  is  presented  in  the  following  section. 


9 . 2  Discussion  of  Results; 

In  the  solution  of  differential  equations  9.1  to  9.5  it  was  assumed 
that  the  incident  flux  is  monochromatic  in  nature.  This  considerably  simplifies 
the  problem  of  selective  absorption  of  radiation  by  the  propellant.  Such  a 
situation  practically  exists  in  ignition  experiments  using  laser  beams.  Since 
the  gas  phase  reactions  are  almost  negligible  up  to  ignition,  there  is  a 
tremendous  build  up  of  reactant  concentration  and  once  the  exothermic  reactions 
are  triggered,  an  enormous  amount  of  energy  is  released  in  the  gas  phase  almost 
instantaneously.  Chemical-diffusional  times  are  extremely  short  compared  to 
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the  overall  ignition  time.  For  example,  in  a  typical  case,  the  chemical  delay 
time  is  on  the  order  of  a  millisecond  whereas  the  total  delay  period  is  on 
the  order  of  200  milliseconds.  Hence  the  ignition  criterion  based  on  a 
certain  surface  temperature  rise  seems  to  be  quite  realistic.  Only  at  extremely 
high  rates  of  heating,  does  the  system  of  gas  phase  equations  really  need  to 
be  considered.  At  high  heating  rates,  the  ignition  behaviour  under  arc  image 
furnace  conditions  should  be  very  much  similar  to  shock  tube  ignition  data. 

At  very  low  heating  rates  (Q  *£3t  5  to  10  cal/cm  sec),  conduction  in  the 
solid  becomes  an  important  factor  and  again  conventional  methods  can  be  employed 
to  predict  the  ignition  behaviour.  All  in  all,  it  appears  that  although  gas 
phase  reactions  play  an  important  role  in  establishing  a  gas  phase  flame,  as  far 
as  the  ignition  behaviour  is  concerned,  gas  phase  processes  can  be  considered 
secondary.  Pressure  seems  to  have  little  effect  on  ignition  behaviour  except 
at  sub-atompsheric  ,or  atmospheric  pressures.  For  example,  reference  to  Figures 
(32),  (33)  and  (34)  reveal  that  above  25  atmospheres,  the  effect  of  any  further 
increase  in  the  gas  pressure  has  virtually  no  effect  on  ignition  delay. 

It  is  interesting  to  note  that  the  exact  numerical  calculations  of  the 
present  investigations  are  in  excellent  agreement  with  the  calculations  performed 
in  reference  (18)  using  approximate  methods.  It  was  also  found  that  at  low  or 
moderately  high  heating  rates  and/or  at  low  pressures,  the  effect  of  initial  gas 
phase  oxidizer  mole  fraction  had  no  effect  cn  the  ignition  behaviours.  However, 
this  may  not  be  the  case  at  extremely  high  heating  rates. 

The  effect  of  absorption  of  radiation  of  the  gas  phase  media  upon  the 
ignition  behaviour  is  very  little.  This  is  understandable  because  radiation 
passess  through  the  gaseous  medium  virtually  unabsorbed,  except  at  very  short 
,.'>ve  lengths.  Even  a  grey  gas  absorbs  very  little  radiation  in  the  visable  and 
infrared  regions  unless  it  passes  through  an  extremely  dense  medium  or  the 
radiation  has  to  be  passed  through  an  extremely  thick  region. 

In  the  present  calculations,  reflection  of  the  radiation  from  the 
propellant  surface  is  not  considered.  Since  the  reflection  of  radiation  from 
the  propellant  surface  depends  entirely  upon  the  surface  roughness,  opacity 
etc.,  the  effect  of  considering  the  reflectivity  ("f)  is  equivalent  to  consider¬ 
ing  an  incident  effective  flux  equal  to  Qo(l-lf)» 

9.3  Conclusions: 


The  present  investigations  has  revealed  that  in  sofar  as  slow  ignition 
is  considered,  like  ignition  in  arc  image  furnaces,  the  simple  theory  based  on 
conduction  and  absorption  of  energy  in  solid  phase,  neglecting  the  gas  phase, 
can  be  used  to  predict  the  ignition  behaviour  of  solid  propellants  at  high  and 
moderate  pressures.  Only  at  very  low  pressures  —  in  fact,  sub-atmospheric 
pressures  —  when  chemical-diffusion  reaction  times  become  very  large,  does  the 
overall  ignition  time  depend  strongly  upon  phase  chemical  kinetics. 

9 . 4  Suggestions  for  Future  Work; 

The  gas  phase  theory  of  ignition  treated  so  far  can  successfully  explain 
several  aspects  of  propellant  ignition  found  during  experiments.  It  is  necessary 
for  successful  application  of  the  present  theory  that  more  precise  experiments 
be  performed.  Without  sensitive  experiments,  it  is  less  meaningful  to  develop  a 
mole  precise  theoretical  model.  As  said  in  previous  chapters,  strong  emission 
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may  be  occurring  close  to  the  establishment  of  steady  state.  To  detect 
light  emission  in  early  stages,  it  is  advisible  to  employ  a  photo-multiplyer 
•rtiich  is  extremely  sensitive. 

Further  improvement  in  the  theory  can  be  brought  about  by  consider¬ 
ing  small  spherical  oxidizer  crystals  embeded  in  fuel  matrix  and  writing  the 
appropriate  differential  equations  and  boundary  conditions.  This  represents 
a  mole  realistic  composite  propellant  mouel. 

It  is  also  of  interest  to  consider  theoretically  the  ignition  of  a 
heterogeneous  propellant  model  by  radiant  energy.  These  combined  with 
sensitive  experimental  should  hopefully  lead  to  a  better  understanding  of 
the  ignition  behaviour  of  solid  propellants. 
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TABLE  1 


Curve  No. 


Description 


1  Long  time  data,  mQ  =  0.05,  0.005 

E/R  =  12000,  e  =  0.01 

2  Long  time  data,  mQ  =  0.05,  0.005 

E/'  =  6000,  e  =  0.01 

3  Short  time  data,  mQ  =  0.005,  E/R  =  6000 

[L*(0)  =  1.56o] 

4  Short  time  data,  mQ  =  0.05,  E/R  =  6000 

(L*(0)  =  1.50o) 


TABLE  2(a) 

CURVE  DESCRIPTION  FOR  FIGURE 


No .  Curve 


Reference 


1  ONERA  6 

20%  PBAA  +  80%  NH, CIO. 

4  4 

2  McAlevy  14 

(Double  base) 

3  McAlvey  6 

(22.5%  Polyester  +  77.5%  NH^XIO^) 

4  Emission  Criterion, 

J*  =  0.1 

5  ’Long  Time  Criterion' 

(e=0.01,  note  the  values  on  the 
curve  are  to  be  multiplied  by  10) 

6  Emission  Criterion 

J*  =  0.01 


v 


TABLE  2(b) 


CURVE  DESCRIPTION  FOR  FIGURE 


No. 

Curve 

Reference 

1 

McAIevy 

(20%  Epoxy  __  80%  NH^Clo^) 

6 

2 

'Light  Emission  Criterion' 

J*  =  0.1,  E/R  =  12000°K 

3 

Kurylko,  L. 

(20%  PBAA  +  80%  NH4C10a 

6 

A 

emission  Criterion 

J*  =  0.01,  E/R  =  6000°K 

TABLE  3 

ASSUMED  VALUES  OF  PHYSICAL  PARAMETERS  USED  IN  THEORETICAL  CALCULATIONS 


Quantity 

Value 

mQ,  g/c.m^-sec  at  400°K 

0.05  and  0.005 

c  ,  cal/g-°K 

0.3 

X  ,  X  ,  cal/cm-sec-°K 
s  g’ 

0.005 

h  ,  cal/gm 

V 

120 

„  2. 

D,  cm  /sec 

0.1 

T  ,  °K 
s,o 

300 

o 

o 

1800 

=  Q2,  cal/mole 

12000 

D  /R,°K 

P 

6000 

TABLE  4 
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FIGURE  1 


SKETCH  OF  NEAR  STEADY  STATE  SPATIAL 
DISTRIBUTION  OF  FUEL,  ,  OXIDIZER, 
rjv  AND  TEMPERATURE,  0,  IN  SECONDARY 
REACTION  ZONE 
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FIGURE  7 


IGNITION  DELAY,  t  ,  msec 


EFFECT  OF  CHOICE  OF  IGNITION  CRITERION  UPON 

IGNITION  DELAY 
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FIGURE  9 
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EFFECT  OF  INITIAL  OXIDIZER  MOLE 
FRACTION  ON  THE  IGNITION  DELAY 
PREDICTED  BY  THE  FEEDBACK  MODEL 
USING  DIFFERENT  IGNITION  CRITERIA 


INITIAL  OXIDIZER  MOLE 
FRACTION  ,  y®x 


FIGIJRF  10 


OPTICAL  DETECTION  SYSTEM 
(  From  reference  I ) 
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PHOTO  CELL 


RELATIVE  THEORETICAL  RESPONSE,  J,  OF  A  LIGHT  EMISSION  DETECTOR 
VERSUS  TIME. 


FIGURE  \?. 


OSCILLOSCOPE  DEFLECTION,  CM  RELATIVE  DETECTOR  RESPONSE,  J 
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THEORETICAL  AND  EXPERIMENTAL  CURVES  OF  THE 
RESPONSE  OF  A  LIGHT  EMISSION  DETECTOR  VERSUS 
TIME  IN  THE  FORMAT  OF  EXPERIMENTAL  DATA. 


(a) 


TIME,  MSEC 
(b) 


FIGURE  13 


EFFECT  CF  INITIAL  OXIDIZER  MOLE  FRACTION  UPON  IGNITION  BEHAVIOUR 
WITH  DIFFERENT  LIGHT  EMISSION  IGNITION  CRITERIA. 
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COMPARISON  OF  THEORETICAL  IGNITION  BEHAVIOR,  USING  VARIOUS 
IGNITION  CRITERIA,  WITH  EXPERIMENTAL  RESULTS. 
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EFFECT  OF  INITIAL  OXIDIZER  MOLE 
FRACTION  ON  THE  THEORETICAL 
IGNITION  DELAY,  FOR  PROPELLANTS 
AND  PURE  FUELS,  USING  VARIOUS 
LIGHT  EMISSION  CRITERIA 
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FIGURE  16 


THEORETICAL  OSCILLOSCOPE  TRACES  FOR  A  HO/iOGENEOlS 

SOLID  PROPELLANT 
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EFFECT  OF  GAS  PHASE  TOTAL 
PRESSURE  ON  IGNITION  DELAY  FOR 
A  HOMOGENEOUS  SOLID  PROPELLANT 


FIGURE  18a 
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EFFECT  OF  GAS  PHASE  TOTAL 

PRESSURE  ON  IGNITION  DELAY  FOR 
A  HOMOGENEOUS  SOLID  PROPELLANT 
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FIGURE  18  b 
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EFFECT  OF  OXYGEN  PARTIAL  PRESSURE  ON 
IGNITION  DELAY  FOR  A  HOMOGENEOUS  SOLID 
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BEHAVIOR  OF  THE  IGNITION  DELAY 
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FIGURE  19 
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Figure  20 


GEOMETRY  OF  A  HETEROGENEOUS 
PROPELLANT 


HATCHED  POSTIONS  REPRESENT  OXIDIZER  PARTICLES 
AND  THE  DOTTED  LINES,  THE  BOUNDARY  OF  FUEL 

SURROUNDING  IT. 


Figure  21 
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igure  24 


THEORETICAL.  "OSCILLOSCOPE"  TRACES  OF 
RELATIVE  LIGHT  EMISSION  VERSUS  TIME 
FOR  HETEROGENOUS  SOLID  PROPELLANT 
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Figure  26 


EFFECT  OF  INITIAL  OXIDIZER  MOLE  FRACTION  ON 
IGNITION  DELAY 

(HETEROGENEOUS  SOLID  PROPELLANT) 
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Figure  27 


RELATIVE  RESPONSE,  Jx5C  RELATIVE  RESPONSE,  Jx50 
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THEORETICAL  OSCILLOSCOPE  TRACES  FOR  A  HETEROGENEOUS 

SOLID  PROPELLANT 
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EFFECT  OF  GAS  PHASE  TOTAL 
PRESSURE  ON  IGNITION  DELAY  FOR 
A  HETEROGENEOUS  SOLID  PROPELLANT 
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FIGURE  29 
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EFFECT  OF  GAS  PHASE  TOTAL 
PRESSURE  ON  IGNITION  DELAY  FOR 
A  HETEROGENEOUS  SOLID  PROPELLANT 
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FIGURE  30 
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FIGURE  31 
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EFFECT  OF  RADIANT  ENERGY 
FLUX  UPON  IGNITION  DELAY 
OF  A  HOMOGENEOUS  PROPELLANT 
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FIGURE  33 


EFFECT  OF  GAS  PHASE  PRESSURE  ON 
THE  IGNITION  DELAY  OF  HOMOGENEOUS 
PROPELLANTS  DURING  RADIANT  IGNITION 
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